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A SEISMIC  RISK  SIMULATION  MODEL  FQR  ARMY 
FACILITIES:  PHASE  ONE,  DEVELOPMENT  OF 
DETERMINISTIC  MODEL 

1 INTRODUCTION 

Background 

Because  the  lateral -force-resisting  systems  of  existing  Army  facil- 
ities were  designed  for  seismic  loads  significantly  less  than  cur- 
rent requirements,  such  facilities  may  provide  inadequate  protection 
from  seismic  loads  defined  in  current  specifications  and  thus  may  pose 
a hazard  to  life  and  property.  To  assess  this  hazard,  existing  facil- 
ities must  be  evaluated  for  seismic  risk.  The  evaluation  must  provide 
the  information  needed  to  determine  what  action  to  take  to  mitigate 
the  potential  hazard. 

Three  major  courses  of  action  may  be  considered:  the  existing 
facility  may  be  (1)  left  unchanged,  with  damage  resulting  from  seismic 
loads  being  repaired  at  some  future  date,  (2)  strengthened  to  meet 
current  specifications,  or  (3)  replaced  with  a facility  meeting  cur- 
rent design  specifications.  Thus,  the  problem  becomes  one  of  pro- 
viding information  that  will  assist  in  making  a rational  decision  re- 
lative to  these  three  courses  of  action. 

Purpose 

The  purpose  of  this  study  is  to  develop  a decision  tool  for  as- 
sessing (1)  the  seismic  hazard  to  Army  facilities  and  (2)  the  cost  of 
mitigation  schemes  for  reducing  the  hazard.  This  report  describes  the 
first  phase  of  the  study,  in  which  a computer  simulation  model  was  devel- 
oped. The  model  simulates-'t'fiO'esf^  bl'~a”facTrity7  predkU  the  damage 
to  be  expected  for  a given  level  of  seismic  activity,  and  develops 
cost  information  for  the  three  decision  alternatives  defined  above. 

Future  phases  of  development  will  contain  refinements  that  will  im- 
prove the  model's  sensitivity  and  provide  additional  decision  informa- 
tion. 
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Approach 

Seismic  Risk  Analysis  Overview 

Underlying  the  process  of  providing  the  required  decision  infor- 
mation is  the  technology  of  seismic  risk  analysis.  Seismic  risk 
analysis  is  a set  of  activities  that  results  in  determination  of  the 
probability  of  occurrence  of  a certain  level  of  seismic  activity  at  a 
particular  place  and  time.  The  risk  analysis  overview  presented  in 
this  section  is  tailored  to  the  U.S.  Army  Corps  of  Engineers  construc- 
tion program  and  is  in  part  taken  from  notes  provided  to  the  Corps  by 
H.  C.  Shah  of  Stanford  University.* 

The  assessment  of  risk  to  an  existing  facility  from  seismic  ac- 
tivity can  be  considered  in  two  parts.  The  "expected  seismic  activity" 
at  the  facility  site  is  determi-ed  first,  followed  by  the  effects  of 
this  expected  seismic  activity  on  the  facility.  The  information  from 
these  two  steps  is  used  in  a decision  methodology  to  provide  guidance 
on  either  accepting  or  reducing  the  risk. 

Determining  the  expected  seismic  activity  at  a site  requires 
that  a site  risk  analysis  be  performed.  In  simplest  terms,  the  site 
risk  analysis  involves  modeling  historical  seismic  data  for  the  region 
to  derive  expressions  for  the  future  expected  seismic  activity.  Shah 
identifies  four  basic  steps  in  determining  the  seismicity  of  a site: 

1.  Determination  of  the  location  of  seismic  sources  by  combining 
seismological  (e.g.,  attenuation)  and  geological  (e.g.,  location  and 
description  of  faults)  information  for  a selected  area  surrounding  the 
site  with  information  on  past  seismic  events  affecting  the  site.  These 
data  provide  basic  qualitative  information  on  the  susceptibility  of  the 
site  to  seismic  activity. 

2.  Description  of  the  seismicity  of  the  sources.  Included  in 
this  step  are  the  models  for  forecasting  magnitude,  intensity, 


H.  C.  Shah,  "Seismic  Risk  Analysis,"  Notes  for  course  given  to  Corps 
of  Engineers,  August  11-15,  1975,  San  Francisco,  CA. 
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duration,  and  peak  ground  acceleration  at  a site  based  on  the  informa- 
tion from  the  first  step.  These  models  are  usually  simple  probabil- 
istic models  that  construct  histograms  for  the  site  based  on  past  data 
and  attenuation  relationships.  A computer  program  that  will  perform 
the  histogram  calculations  based  on  a simple  Poisson  pr  v.-ilistic 
model  is  available.^ 

3.  Careful  study  of  the  microcharacteristics  of  the  site,  in- 
cluding characteristics  of  the  soil-structure-foundation  interaction. 
Although  this  interaction  is  not  well  understood,  it  is  often  of  cri- 
tical importance. 

4.  Assembling  data  to  provide  information  on  peak  ground  acceler- 
ation at  the  site,  and  duration  and  number  of  occurrences  of  a parti- 
cular level  of  acceleration.  For  any  computerized  seismic  simulation 
model,  knowledge  of  the  peak  ground  acceleration  is  essential.  Knowl- 
edge of  the  number  of  occurrences  is  essential  for  the  decision  method- 
ology, if  other  than  a deterministic  approach  to  seismic  risk  assess- 
ment is  taken. 

The  second  part  of  the  risk  assessment  problem  considers  the  ef- 
fects of  seismic  activity  on  the  facility.  This  study  takes  an 
analytical  approach  to  prediction  of  damage— perhaps  the  most  straight- 
forward approach  for  those  inexperienced  in  seismic  hazard  evaluation. 
It  is  based  on  elastic  code  design  of  the  facility;^  i.e.,  nonlinear 
effects  are  neglected.  Structural  behavior  under  code-defined  loading 
is  used  to  predict  the  damage  to  be  expected  from  a given  level  of 
seismic  activity.  The  analysis  model  is  a simple  story  stiffness  type. 
The  total  simulation  model  goes  beyond  damage  prediction  and  provides 
guidance  for  upgrading  needed  to  bring  the  facility  within  acceptable 
risk  limits. 


" H.  C.  Shah. 

^ Seismic  Design  for  Buildings,  TM  5-809-10  (Department  of  the  Army, 
April  1973);  and  Recommended  Lateral  Force  Requirements  and  Com- 
mentary (Seismology  Committee,  Structural  Engineers  Association  of 
California  [SEAOC],  1968). 
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simulation  Model  Requirements 


; 


Model  simulation  is  a commonly  used  technique  in  decision  method- 
ology. Construction  of  a simulation  model  requires  that  the  assumptions 
be  clearly  presented  and  that  all  factors  of  significance  be  accounted 
for.  Once  the  model  is  constructed,  it  must  be  capable  of  transforming 
the  facility  configuration  into  information  for  use  in  deciding  whether 
to  repair  the  facility  after  an  earthquake,  to  strengthen  the  facility, 
or  to  replace  it.  Facility  configuration  input  must  be  minimal;  only 
information  pertinent  to  the  decision  should  be  required.  The  model 
must  provide  information  within  a short  time  frame  and  be  flexible 
enough  to  be  reused  at  minimal  cost  whenever  new  requirements  are  pro- 
posed. All  facilities  must  be  compared  on  the  same  basis,  and  differ- 
ences in  decision  information  must  result  from  differences  in  the 
facility  configurations. 

The  model  described  in  this  report  was  developed  to  meet  the  above 
requirements. 
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2 SEISMIC  RISK  SIMULATION  MODEL 

t 

The  seismic  risk  simulation  model  translates  the  configuration  of 
a given  facility  into  information  pertinent  to  assessing  the  risk  to 
the  facility  from  seismic  loads.  The  model  analyzes  four-  to  ten- 
story-high  reinforced  concrete  structures  having  a reasonably  sym- 
metric configuration.  The  model  consists  of  four  parts,  which  can  be 
referrecf^to  as  input,  design,  analysis,  and  cost  tabulation.  Appen- 
dix A describes  the  model  in  detail,  while  this  chapter  defines  its 
scope  and  limitations. 

Input 

Input  to  the  model  includes: 

1.  Number  of  stories 

2.  Story  height  (assumed  same  for  all  floors) 

3.  Story  weight 

4.  Number  of  bays  in  each  direction 

5.  Ratios  of  moments  of  inertia  of  shear  wall  to  frame  columns 
over  the  building  height 

6.  Overturning  moment  factor 

7.  Floor  area  per  column 

8.  Total  construction  cost  per  square  foot  for  structural  ele- 
ments 

9.  Total  construction  cost  per  square  foot  for  nonstructural 
elements 

10.  Structure  loads  in  terms  of  code  variables 

11.  Convergence  option 

12.  Calibration  constants. 

For  an  existing  facility,  most  of  the  information  in  items  1 
through  10  is  readily  available.  Except  for  story  height  and  struc- 
tural and  nonstructural  construction  costs,  all  variables  are  story- 
dependent,  allowing  for  flexibility  in  the  specification  of  the 
structural  configuration.  A static  analysis  is  performed  on  the 
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elastic  model,  and  loads  on  the  structure  are  defined  in  terms  of  the 
code  used  to  design  the  facility. 

Chapter  3 describes  the  procedure  for  determining  the  calibration 
constants,  which  correlate  cost  information  with  the  variables  related 
to  the  structure's  behavior. 

Design 

The  model  simulates  the  design  of  the  structure;  that  is,  the 
program  sizes  the  column  elements  and  shear  wall  elements  so  that  their 
resistance  meets  or  exceeds  the  seismic  code  requirements.  The  simu- 
lated design,  which  only  approximates  a detailed  design,  is  necessary 
to  determine  parameters  to  strengthen  and  replace  the  facility.  Since 
a priori  rationales  for  proceeding  from  the  existing  facility  to  a 
strengthened  or  suitable  replacement  facility  do  not  exist,  simulation 
of  the  facility  design  is  accomplished  using  the  following  algorithm. 

An  initial  frame  column  depth  is  assumed.  From  this,  the  stiffness 
matrix  for  the  frame  is  generated.  If  shear  walls  are  to  be  included 
in  the  configuration,  the  stiffness  matrix  for  the  shear  wall  elements 
is  also  determined.  The  frame  and  shear  wall  stiffness  matrices  are 
added,  the  matrix  inverted,  and  story  displacements  determined  along 
with  the  moments  and  forces  in  the  frame  members. 

If  the  member  moments  and  forces  exceed  the  capacity  of  the  ele- 
ment, the  size  of  the  element  (and  thus  its  capacity)  is  increased; 
if  the  moments  and  forces  within  the  frame  element  are  less  than  the 
capacity  of  the  element,  the  size  of  the  element  may  be  decreased, 
depending  on  the  convergence  option.  This  process  is  repeated  until 
the  member  forces  are  less  than  the  member  capacity.  Convergence  for 
a frame  structure  is  obtained  within  a few  iterations.  If  shear  walls 
are  included  and  if  convergence  is  not  obtained  within  a specified 
number  of  iterations,  the  design  is  assumed  not  to  be  adequate. 

Factors  entering  into  the  design  include  dead  load,  live  load, 
seismic  load,  and  overturning  moment.  To  simplify  the  design  proce- 
dure, the  floor  diaphragms  and  their  accompanying  beams  and  girders 
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are  assumed  infinitely  rigid,  and  columns  are  assumed  to  have  fixed 
ends.  ’ 


Analysis 

The  analysis  portion  of  the  model  determines  the  member  forces 
from  the  computed  story  displacements.  The  member  forces  are  based 
on  second-order  displacement  difference  expressions. 

Cost  Tabulations 

The  model  also  provides  cost  computations.  The  following  para- 
graphs provide  an  overview  of  the  three  cost  algorithms  in  the  model. 
Whenever  cost  is  referred  to,  it  is  assumed  that  the  cost  is  the  original 
construction  cost  for  the  structure. 

The  damage/repair  algorithm  is  based  on  calculation  of  a damage 
factor.  From  Blume's  work,**  the  reserve  energy  of  a structure  is 
equated  to  the  area  under  a hypothetical  elastoplastic  load  deflection. 
This  curve  provides  an  expression  for  ductility  in  terms  of  the  demand 
placed  on  the  structure  relative  to  its  capacity  to  resist  that  de- 
mand. From  this  expression,  a damage  factor  which  is  a measure  of 
the  reserve  capacity  of  a structure  may  be  computed.  A portion  of 
the  structural  and  nonstructural  damage  resulting  from  seismic  loading 
greater  than  the  facility  design  level  is  assumed  proportional  to  the 
damage  factor. 

It  is  generally  accepted  that  nonstructural  damage  is  related  to 
factors  other  than  the  overall  damage  factor.  One  such  factor  is  re- 
lative story  displacement;  i.e.,  increased  relative  story  displacement 
correlates  positively  with  increased  nonstructural  damage.  An  addi- 
tional expression  for  nonstructural  damage  is  therefore  included  in 
the  model. 

Thus,  the  total  cost  to  repair  the  structure  is  the  sum  of  the 
general  structural  and  nonstructural  damage  factor  for  the  overall 
structure  and  the  nonstructural  damage  factors  for  the  individual 

Ef feats  Prediction  Guidelines  for  Struotures  Subjected  to  Ground 
f Motion,  JAB-99-115,  UC-11  (URS/John  A.  Blume  and  Associates,  Engi- 

neers,  July  1975). 

h 
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stories.  Although  it  is  recognized  that  other  factors,  such  as  struc- 
ture torsion  and  foundation  characteristics,  inay  also  contribute  to 
substantial  structural  and  nonstructural  dama^jH,  these  factors  are 
not  accounted  for  in  this  initial  model. 

The  strengthening  algorithm  consider*"  ‘he  d.fference  in  stiff- 
nesses between  the  designs  simulated  for  tne  current  requirements  and 
the  requirements  under  which  the  facility  was  designed.  As  a criter- 
ion for  strengthening,  the  story  dispel acemerts  of  the  strengthened 
design  must  be  less  than  or  equal  to  those  of  the  initial  design.  When 
shear  walls  are  used  to  strengthen  the  structure,  the  diagonal  terms 
of  the  stiffness  matrix  are  representative  of  relative  story  stiff- 
nesses and  are  used  to  compute  the  strengthening  costs.  For  nonshear 
wall  cases,  the  stiffness  matrix  to  strengthen  the  structure  is  gener- 
ated initially  by  taking  10  percent  of  the  initial  frame  design  stiff- 
ness matrix.  The  total  sfengthening  cost  for  both  the  structural  and 
nonstructural  elements  is  a function  of  the  difference  between  the 
stiffness  of  the  existing  facility  and  that  of  the  simulated  design  of 
the  strengthened  facility. 

The  replacement  algorithm  measures  the  increase  m material  needed 
in  the  simulated  design  for  higher  requirements.  The  column  area  and 
the  percentages  of  shear  wall  steel  difference  provide  a measure  of 
the  increase  in  material  costs  due  to  replacement.  The  structural 
costs  are  considered  proportional  to  these  ratios.  The  nonstructural 
costs,  which  are  not  likely  to  increase  under  replacement,  are  con- 
sidered proportional  to  a calibration  factor  times  the  original  non- 
structural costs. 
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3 MODEL  CALIBRATION 

Model  calibration  is  the  process  of  determining  unknown  model  con 
stants  which  insure  that  the  information  produced  by  the  model  agrees 
with  available  data.  The  process  is  similar  to  curve  fitting.  The 
constants  for  the  model  were  determined  such  that  damage/repair, 
strengthening,  and  replacement  costs  can  be  determined  for  four- 
through  ten-story  reinforced  concrete  facilities  with  a variety  of 
configurations  for  a variety  of  levels  of  seismic  load. 

Table  1 gives  the  matrix  of  building  configurations  under  consid- 
eration. The  building  configurations  are  defined  in  terms  of  three 
parameters— structural  system,  dimensional  ratio,  and  ratio  of  shear 
wall  to  frame  moments  of  inertia.  In  this  study,  the  shear  wall 
moment  of  inertia  is  held  constant  over  the  story  levels,  with  the 
ratio  of  wall  moment  of  inertia  to  frame  moment  of  inertia  being  com- 
puted at  the  first  story.  An  actual  facility  is  transformed  into  a 
story  stiffness  model  by  specifying  the  ratio  of  wall  to  frame  stiff- 
ness for  the  first  story,  and  the  relative  frame  stiffness  at  stories 
above  the  first.  The  structural  systems  range  from  a frame  (K  equals 
0.67),  to  a light  shear  wall  (K  equals  0.80),  to  a massive  shear  wall 
(K  equals  1.33);  the  ratio  of  the  shear  wall  moment  of  inertia  to  the 
frame  moment  of  inertia  varies  from  0 to  5,000.  The  dimensional  ratio 
pertains  to  the  overturning  moment  factor  in  the  model.  The  larger 
the  dimensional  ratio,  the  wider  the  building  in  the  direction  of  the 
load.  The  floor  and  roof  weights  used  in  the  simulation  are  typical 
of  industrial  buildings.  For  computation  of  cost,  14  configurations 
and  a continuum  of  loads  ranging  from  1968  Structural  Engineers  Associ 
ation  of  California  (SEAOC)  Zone  3 requirements  to  a 1.0  g response 
spectrum  ® were  considered. 

In  the  calibration  of  the  model,  data  pertinent  to  the  building 
configurations  in  Table  1 and  the  results  of  initial  run^  "or 

W.  K.  Stockdale,  Modal  Analysis  Methods  in  Seismu  Design  for 
Buildings,  Technical  Report  N-132/ADA012732  (U.S,  Army  Construc- 
tion Engineering  Research  Laboratory  [CERL],  June  1975). 
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Table  1 

Model  Reinforced  Concrete  Building  Facility  Configurations 


[figuration 

Number 

Str^jotural 

oj'stem 

Dimensional 

Ratio 

Ratio  of  Shear 
Wall  Moments  of 
Inertia  to  Frame 
Moment  of  Inertia 
(IW/IF) 

1 

K=0.67 

1.0 

0 

2 

K=0.67 

1.6 

0 

3 

K=0.80 

1.0 

100 

4 

K=0.80 

1.0 

200 

5 

K=0,80 

1.0 

500 

6 

K=0.80 

1.5 

100 

7 

K=0.80 

1.6 

200 

8 

K=0.80 

1.6 

500 

9 

K=1.33 

1.0 

1000 

10 

K=1.33 

1.0 

2000 

11 

K=l,33 

1.0 

5000 

12 

K=1.33 

1.6 

1000 

13 

K=1.33 

1.6 

2000 

14 

K=1.33 

1.6 

5000 
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performance  of  the  model  with  arbitrary  calibration  constants  were 
used  to  determine  a f’nal  set  of  model  calibration  constants  applicable 
to  all  configurations  for  the  above  loading  range.  Appendix  B presents 
the  details  of  the  calibration  process  along  with  the  supporting  data. 

A difficulty  experienced  in  calibration  of  tne  model  was  that 
the  data  available  for  calibration  are  inadequate  for  the  degree  of 
detail  incorporated  into  the  model.  The  damage  study  of  the  San 
Fernando  earthquake  performed  at  the  Massachusetts  Institute  of  Tech- 
nology (MIT),®  which  provided  much  of  the  supporting  data,  does  not 
consider  the  structural  system  type.  Strengthening  data  are  sparse 
and  keyed  to  critical  facilities.  Replacement  data  are  generally 
available  only  for  moderately  low  seismic  levels;  extrapolation  to 
the  levels  considered  in  this  report  may  be  misleading. 


T R.  V.  Whitman  et  a1..  Damage  Statistics  for  High-Rise  Buildings 
in  the  Vicinity  of  the  San  Fernando  Earthquake^  R73-24,  Report  No. 
7,  Structures  Publication  363  (Department  of  Civil  Engineering, 
Massachusetts  Institute  of  Technology  [MIT],  April  1973). 
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l^  EXAMPLE  MODEL  OUTPUT 

This  chapter  presents  an  example  of  the  initial  information  ob- 
tained from  the  model  for  a variety  of  structural  configurations. 

For  purposes  of  illustration,  it  is  assumed  that  a number  of  rein- 
forced concrete  building  configurations  (four  to  ten  stories)  are  to 
be  assessed  for  seismic  hazard  and  determination  of  the  cost  of  the 
hazard  mitigation  schemes  presented  in  Chapter  1.  It  is  also  assumed 
that  the  facilities  safsfy  1968  SEAOC  Zone  3 requirements,  but  now 
must  meet  higher  requirements  (up  to  and  including  a 1.0  g response 
spectrum  acceleration). 

To  obtain  cost  information  on  a particular  building  configura- 
tion, the  configuration  parameters  required  to  create  a story  stiffness 
model  can  be  provided  and  the  simulation  model  run.  This  is,  however, 
time  consuming,  particularly  if  several  hundred  buildings  are  under 
consideration.  To  obtain  information  on  a particular  building  con- 
figuration easily  without  rerunning  the  model  for  each  particular  con- 
figuration, representative  building  configurations  are  selected,  the 
cost  information  determined,  and  interpolations  of  the  cost  curves  ac- 
cording to  configurations  employed.  Thus  one  set  of  cost  curves 
determined  for  a set  of  representative  building  configurations  provides 
information  on  many  particular  building  configurations. 

Of  the  42  configurations  (four,  seven,  and  ten  stories)  in  Table  1, 
21  were  selected  for  determination  of  cost  information.  The  dimension- 
al ratio  was  assumed  to  be  1.6,  and  the  ratio  of  shear  wall  to  ^rame 
stiffness  assumed  to  be  0,  100,  200,  500,  1000,  2000,  or  5000  for  four-, 
seven-,  and  ten-story  configurations.  This  set  of  21  structural  con- 
figurations provides  the  building  configuration  data  space.  Costs  to 
repair,  strengthen,  and  replace  were  derived  for  these  configurations 
for  seismic  loads  ranging  from  about  0.3u  g to  1.0  g response  spectra. 
Limited  extrapolation  will  allow  for  somewhat  higher  and  lower  levels. 

Figures  1 through  9 present  the  configuration-load-cost  data  for 
15  of  the  21  configurations  with  damage/repair,  strengthening,  and 
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replacement  percentage  costs  of  the  original  facility  tabulated  sepa- 
rately. Figures  10  through  15  present  the  percentage  costs  of  the  six 
remaining  facilities  with  damage/ repair  (D/R),  strengthening  (S),  and 
replacement  (R)  costs  on  one  graph.  The  relative  magnitudes  of  costs 
of  the  decision  alternatives  can  be  easily  ascertained  for  each  facil- 
ity. The  facility  numbers  used  in  the  figures  refer  to  the  configu- 
rations listed  in  Table  1. 

The  model  cost  information  was  analyzed  to  facilitate  understand- 
ing of  the  cost  estimates  derived  from  the  model  and  of  the  costing  al- 
gorithms detailed  in  Appendix  A,  as  well  as  to  provide  information  for 
improving  the  cost  algorithms  in  subsequent  refinements  of  the  simula- 
tion model.  This  analysis  led  to  the  following  observations  on  the 
damage/ repair,  strengthening,  and  replacement  costs  for  the  four-, 
seven-,  and  ten-story  structures: 

1.  Damage/repair  costs 

a.  Damage/repair  costs  increase  with  demands  on  the  struc- 
ture for  all  structure  configurations. 

b.  The  model  is  not  as  sensitive  to  damage/repair  as  would 

be  desirable;  damage/ repair  costs  are  near  80  percent  for  all  structure 
configurations.  The  reason  for  this  is  that  the  damage  factor  for 
structural  damage  (DF^)  is  1.0  at  comparatively  low  levels  of  demand 
because  of  the  moderately  low  value  of  ultimate  ductility  (p^  equals 
5.0).  For  a 1.0  g response  spectrum  demand  for  a facility  designed  for 
1968  SEAOC  Zone  3 requirements,  the  damage  should  probably  be  100  per- 
cent, whereas  for  a 0.25  g response  spectrum  demand  on  the  same  facil- 
ity, the  damage  could  be  less  than  one-half  of  this  value. 

c.  The  story  displacement  damage  factor  (OF^^^)  helps  in- 
crease the  sensitivity  of  the  model  with  respect  to  damage/repair 
costs.  In  refinements  of  this  model-,  damage  may  be  defined  only  in 
terms  of  relative  story  displacement. 

d.  Damage/repair  costs  generally  decrease  with  increasing 
ratios  of  wall  stiffness  to  frame  stiffness.  This  seems  reasonable  be- 
cause of  the  corresponding  decrease  in  relative  story  displacement. 
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e.  Two  important  variables  aid  in  understanding  the  cost 
informat ion- -mean  cost  and  standard  deviation  of  cost.  The  first  var- 
iable provides  a measure  of  the  relative  magnitude  of  the  cost  data, 
whereas  the  second  variable  determines  the  dispersion  of  the  cost  data. 
Table  2 presents  this  information  for  all  21  configurations.  In  gen- 
eral, for  damage/repair  costs,  the  standard  deviation  (the  dispersion) 
is  least  for  massive  shear  wall  configurations  (K  equals  1.33),  some- 
what greater  for  frame  configurations  (K  equals  0.67),  and  largest  for 
light  shear  wall  configurations  (K  equals  0.80).  In  engineering  terms, 
this  indicates  that  damage/repair  cost  predictions  based  on  this  model 
are  least  uncertain  for  massive  shear  wall  configurations,  more  un- 
certain for  frame  configurations,  and  quite  uncertain  for  light  shear 
wall  configurations.  Here  again,  the  relative  story  displacement  dam- 
age factor  is  the  principal  quantity  defining  the  dispersion  or 
uncertainty.  The  mean  of  the  cost  information  with  respect  to  two  con- 
figuration sizes— large  and  small—for  all  story  cases  for  a structural 
system  with  K equal  to  0.67  indicates  that  a small  configuration  will 
provide  higher  damage/repair  costs.  For  the  four-story  configurations 
where  K equals  0.80  or  1.33,  the  smaller  configuration  in  general  has 

a lower  damage/ repair  cost.  For  the  ten-story  configuration  where  K 
equals  0.80  or  1.33,  the  smaller  configuration  has  the  higher  damage/ 
repair  costs.  These  results  are  again  attributable  to  the  relative 
story  displacement  damage  factor. 

f.  Based  on  this  model,  massive  shear  wall  structures  have 
the  lowest  damage/repair  costs  while  light  shear  wall  structures  have 

t 

the  highest  damage/ repair  costs. 

2.  Strengthening  costs 

a.  Strengthening  costs  increase  almost  linearly  with  the 
demands  on  the  structure  for  all  structural  configurations. 

b.  Strengthening  costs  are  widely  dispersed. 

c.  The  strengthening  algorithm  based  on  relative  stiffness 
is  very  sensitive  to  structural  system  type  and  to  load  magnitude. 

For  the  small  amount  of  strengthening  data  available,  the  model  is 
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Table  2 

Mean  and  Standard  Deviation  of  Costs  for  Configurations 
Over  Decision  Alternatives 


Configuration 

Damage/Repair  Strengthening 

Replacements 

Number 

m 

£ 

m 

£ 

m 

Four  Story 

2 

84.0 

1.9 

348.4 

163.4 

118.5 

9.6 

6 

86.6 

3.0 

118.6 

34.8 

119.0 

8.7 

7 

84.9 

2.2 

87.4 

25.2 

113.0 

6.8 

8 

82.1 

1.0 

43.6 

14.3 

106.1 

4.5 

12 

81.1 

.3 

19.9 

14.2 

107.4 

7.6 

13 

80.6 

.1 

7.1 

8.8 

104.2 

6.1 

14 

80.2 

.1 

3.9 

7.9 

104.1 

6.4 

Seven  Story 

2 

82.7 

1.3 

394.7 

187.2 

127.8 

14.4 

6 

89.7 

4.5 

240.8 

107.7 

124.4 

14.2 

7 

89.0 

4.2 

167.7 

60.9 

117.3 

9.3 

8 

86.9 

3.2 

99.2 

29.3 

110.2 

5.9 

12 

83.2 

.9 

52.2 

17.1 

110.4 

7.2 

13 

81.7 

.5 

50.8 

17.0 

110.5 

7.6 

14 

80.7 

.2 

50.1 

16.8 

110.6 

7.6 

Ten  Story 

2 

82.0 

1.0 

412.5 

197.5 

123.9 

12.4 

6 

91.4 

5.0 

409.2 

199.5 

119.0 

11.2 

7 

91.0 

4.9 

300.6 

144.7 

115.8 

9.5 

8 

88.6 

4.1 

147.7 

61.5 

108,9 

5.5 

12 

83.6 

1.0 

65.0 

22.7 

108.6 

5.5 

13 

82.0 

.6 

62.0 

21.7 

108.8 

5.6 

14 

80.8 

.2 

60.0 

21.0 

108.9 

5.8 
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adequate.  The  model  is  calibrated  on  the  basis  of  a massive  shear  wall 
structure  of  eight,  ten,  and  11  stories;  strengthening  costs  increase 
substantially  for  both  the  light  shear  wall  structure  (K  equals  0.80) 
and  the  frame  structure  (K  equals  0.67).  Strengthening  costs  also  in- 
crease fairly  substantially  as  the  number  of  stories  of  the  building 
increases.  Both  of  these  trends  in  strengthening  cost  data  are  rea- 
sonable. Si*'..e  strengthening  cost  data  are  widely  dispersed,  standard 
deviations  are  often  over  one-third  of  the  mean  value  and  may  approach 
one-half  of  the  value.  This  is  not  the  case  for  either  damage/repair 
costs  or  replacement  costs.  Greater  variability  is  demonstrated  in 
frame  and  light  shear  wall  structures  than  in  massive  shear  wall  struc- 
tures in  general.  The  analysis  of  the  strengthening  algorithm  is 
complicated  by  the  simulated  design  procedure.  At  low  levels,  moment 
appears  to  control  the  column  design  whereas  at  high  levels,  shear  ap- 
pears to  control  the  design. 

d.  Based  on  this  analysis,  massive  shear  wall  structures 
require  the  least  strengthening  while  frame  structures  require  sub- 
stantially more,  with  light  shear  wall  structures  between  the  two. 

3.  Replacement  costs 

a.  Replacement  costs  increase  nearly  linearly  with  the  level 
of  demand  on  the  structure  for  all  structural  configurations. 

b.  Replacement  costs  are  less  dispersed  than  strengthening 
costs,  but  more  dispersed  than  damage/repair  costs.  In  general,  re- 
placement cost  dispersion  decreases  as  tiie  ratio  of  shear  wall  to 
frame  stiffness  increases. 

c.  Replacement  costs  for  light  shear  wall  and  frame  struc- 
tures increase  over  those  for  massive  shear  wall  structures;  that  is, 
replacing  a frame  structure  with  a frame  structure  costs  substantially 
more  than  replacing  a massive  shear  wall  structure  with  a massive 
shear  wall  structure.  Replacement  costs  based  on  increased  material 
usage  are  accurate,  and  the  simulated  design  may  provioe  more  ac- 
curate  information  here  than  in  the  strengthening  and  uamage/repair 

. i 

algorithms.  Replacement  costs  are  slightly  higher  for  seven-story 
facilities  than  for  ten-story  facilities  (a  reflection  on  the  avail- 
able calibration  data). 
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5 decision  methodology 


The  computer  simulation  model  is  a decision  tool  which  uses  de- 
sign simulation  to  translate  a structure's  behavior  into  variables 
that  correlate  with  costs.  Associated  with  cost  information  for  re- 
pair, strengthening,  and  replacement  of  a given  facility  is  a ration- 
ale for  choosing  one  of  the  three  following  decision  alternatives: 

1.  Leave  the  facility  unchanged  and  anticipate  repair  of  the 
structure  in  the  case  of  damage  due  to  seismic  activity  (inherently 
one  assumes  that  hazard  to  life— a variable  not  quantified  in  the 
model— is  small  and  acceptable) 

2.  Strengthen  the  facility  by  addition  of  lateral-force-resisting 
elements  to  present  code  levels 

3.  Replace  the  facility  with  a facility  designed  for  present 
code  levels. 

The  rationale  for  making  the  decision  may  be  simple  (minimum 
cost)  or  complex  (utility  theory).  It  is  important  to  realize  that 
the  decision  rationale  is  not  built  into  the  simulation  model,  but  is 
an  adjunct.  This  chapter  describes  two  possible  decision  rationales. 

The  first  decision  rationale  is  termed  the  Current  Minimum  Cost 
Rationale  (CMCR).  This  decision  rationale  states  "select  the  alter- 
native that  will  result  in  the  minimum  cost  based  upon  the  cost  infor- 
mation from  the  computer  model."  This  rationale  would  be  the  most 
logical  if  uncertainty  is  not  considered,  i.e. , if  the  problem  is  de- 
terministic. 

The  second  decision  rationale,  termed  the  Expected  Monetary  Cost 
Rationale  (EMCR),  is  an  extension  of  the  first,  in  that  it  introduces 
a level  of  uncertainty  quantification  into  the  decision  procedure. 

The  EMCR  states  "select  the  alternative  that  will  result  in  the  mini- 
mum expected  monetary  cost  based  upon  the  computer  simulation  model 
and  the  site  dependent  seismic  history."  The  expected  monetary  cost 

for  an  event  E,  is  determined  by  multiplying  the  cost  associated  with 
1 

the  event  (C  ) times  the  probability  of  its  occurrence  (P(E)).  The 
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result  is  interpreted  as  the  amount  of  money  the  event  may  be  expected 
to  cost  over  a long  period  of  time.  If  one  considers  all  possible 
events  in  a mutually  exclusive  way  (i.e.,  no  two  events  have  anything 
in  common  with  each  other),  the  chance  or  probability  of  an  event  mea- 
sures the  likelihood  of  the  event  relative  to  a period  of  time.  Thus, 
this  more  sophisticated  rationale  extends  the  cost  estimates  into  the 
future  ;based  on  data  from  the  past.  There  are  a number  of  ways  of 
using  data  based  on  seismic  activity  in  a particular  area  to  construct 
a "seismic  histogram"  defining  the  relative  frequencies  of  the  largest 
expected  seismic  event  in  a given  number  of  years.  Shah’ -has  developed 
a computer  program  which  uses  earthquake  data  in  terms  of  coordinates 
of  the  epicenter  and  magnitude  in  an  appropriate  attentuation  equation 
to  arrive  at  a seismic  histog=*am  for  a particular  site. 

The  model  is  easily  used  in  conjunction  with  this  second  decision 
rationale.  Once  the  seismic  histogram  is  obtained,  the  probability 
of  a certain  Modified  Mercalli  Intensity  (MMI)  level  earthquake  is 
determined.  The  MMI  can  then  be  correlated  with  peak  ground  acceleration 
and,  based  on  this  expected  peak  ground  acceleration,  the  response  spectrum 
approach  can  be  used  to  determine  code  input  requirements  to  the  simula- 
tion model.  For  a particular  structural  configuration,  the  simulation 
model  with  the  given  series  of  design  input  requirements  will  then 
provide  a series  of  costs  to  repair,  strengthen,  or  replace. 

There  are  a host  of  other  decision  rationales.  However,  these 
two  are  perhaps  the  simplest  and  provide  a starting  place.  Applica- 
tion of  both  these  rach  .ales  for  the  cost  data  provided  in  Chapter  4 
is  straightforward  once  the  site  risk  has  been  determined.  The  com- 
pleted simulation  model  will  provide  the  necessary  input  for  both 
these  rationales  in  addition  to  considering  other  decision  rationales. 
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6 CONCLUSIONS  AND  RECOMMENDATIONS 
Conclusions 

This  report  has  described  a computer  simulation  model  developed 
in  the  first  phase  of  a study  designed  to  provide  a decision  tool  for 
assessing  the  seismic  hazard  to  Army  facilities  and  the  cost  of  miti- 
gation schemes  for  reducing  the  hazard.  Observation  of  the  simula- 
tion model's  behavior  and  the  supporting  data  has  led  to  the  following 
conclusions: 

1.  Cost  information  can  be  correlated  with  structural  behavior 
by  selecting  suitable  variables  related  to  the  structural  behavior. 

In  the  model  design,  base  shear,  relative  story  displacement,  building 
stiffness,  and  the  dimensions  of  building  elements  are  the  variables 
considered.  A calibration  procedure  correlates  cost  information  with 
these  variables. 

2.  The  simulated  design  concept  appears  to  provide  the  best 
means  of  comparing  the  behavior  of  structures  at  different  load  re- 
quirements on  the  same  basis.  This  approach  was  necessary  for  estimat- 
ing strengthening  and  replacement  costs. 

3.  For  the  degree  of  detail  incorporated  into  the  model,  the 
data  available  for  calibration  are  inadequate.  The  MIT  damage  study 
on  the  San  Fernando  earthquake  does  not  consider  the  structural  sys- 
tem type.  The  strengthening  data  are  sparse  and  keyed  to  critical 
facilities.  The  replacement  data  are  for  moderately  low  seismic  levels 
and  extrapolation  to  the  levels  considered  in  this  report  may  be  mis- 
leading. 

4.  The  damage/repair  algorithm  selected  is  insensitive  at  high 
levels  of  seismic  force  because  of  the  comparatively  low  value  of  ul- 
timate structural  ductility.  Damage/repair  costs  increase  with  level 
of  seismic  activity  for  a given  building  configuration.  Damage/repair 
costs  are  least  for  massive  shear  wall  configurations  and  highest  for 
light  shear  wall  configurations. 


,'K 

5. '  Strengthening  costs  increase  almost  linearly  with  seismic 
level.  Strengthening  costs  decrease  with  an  increase  in  ratio  of 
shear  wall  to  frame  stiffness.  Strengthening  costs  also  increase  with 
the  building  story  configuration.  ~ 

6.  Replacement  costs  increase  nearly  linearly  with  seismic  level. 
Replacement  costs  in  general  are  highest  for  frame  building  systems 
and  lowest  for  massive  shear  wall  building  systems.  Replacement  costs 
decrease  with  increasing  ratio  of  shear  wall  to  frame  stiffness. 

Recommendations 

The  following  recommendations  are  made  based  on  the  results  of 
the  first  phase  of  this  study: 

1.  The  simulated  design  procedure  in  the  model  should  be  aug- 
mented with  a more  flexible  interactive  design  procedure  in  the  case 
of  strengthening  and  replacement  algorithms,  if  analytical  prediction 
of  strengthening  and  replacement  is  to  be  a viable  approach. 

2.  A site  risk  model  should  be  incorporated  into  the  overall  model. 

3.  The  decision  methodology  should  be  examined  and  implemented. 

4.  Additional  data  for  calibration  of  the  model  should  be  ob- 
tained either  through  disaster  studies  or  statistical  simulation. 

5.  The  simulation  model  should  be  used  as  a true  simulator  for 
developing  repair/damage,  strengthening,  and  replacement  cost  data 
for  many  structural  configurations.  This  is  a realistic  approach  dic- 
tated by  the  lack  of  data  and  the  resources  involved  in  model  calibra- 
tion for  many  particular  cases. 


43 


APPENDIX  A: 


COMPUTER  MODEL 

This  appendix  describes  the  computer  model  in  detail,  including 
the  structural  design/analysis  algorithms  and  the  cost  computation 
models.  The  first  portion  of  the  appendix  presents  the  derivation  of 
the  equations  for  the  elastic  structural  model.  The  second  portion 
describes  the  design  algorithm,  and  the  third  part  contains  a detailed 
description  of  the  cost  computation  models  for  facility  repair, 
strengthening,  and  replacement. 

Derivation  of  the  Structural  Model 
Equations  Governing  Model  Behavior^ 

The  following  assumptions  were  made  in  deriving  the  structural 
model ; 

1.  All  lateral  force  is  resisted  by  a frame  system  of  a combina- 
tion frame-shear  wall  system. 

2.  The  model  foundation  is  rigid. 

3.  The  floor  diaphragms  are  rigid  with  girders  of  infinite 
stiffness. 

4.  Lateral  forces  are  distributed  according  to  the  SEAOC  pro- 
visions. 

5.  Shear  walls  act  as  short  deep  beams  with  stiffness  components 
in  bending  only. 

6.  Structure  overturning  moment  and  dead/live  load  are  accounted 
for  in  the  analysis  and  assumed  to  increase  the  moments  on  frame  and 
shear  wall  elements. 

7.  Sizing  of  elements  is  based  on  column  effective  depths  and 
the  percentage  of  steel  contained  in  the  shear  wall. 


^ J.  C.  Shang,  Design  of  Combined  Frames  and  Shear  Halls,  Advanced 
Engineering  Bulletin  14  (Portland  Cement  Association,  1973). 
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S\  Elongation  and  contraction  of  the  shear  wall  and  column  ele- 
itien^ts  is  neglected. 

The  derivation  is  based  on  a finite  difference  model  of  a struc- 
ture of  N stories  (with  N being  greater  than  or  equal  to  4)  and  three 
bays  in  the  direction  of  the  applied  loads.  This  is  a standard  model  for 
finite  difference  analysis  of  building  structural  models. 

Figure  A1  shows* the-frame  system  configuration  used  in  determin- 
ing the  lateral  forces  distribution.  Equilibrium  of  joint  i requires 

■f  f 

that  the  external  forces  (F'.)  be  balanced  by  the  internal  shears  (V, 

and 


F^  = 

^i  M M'+l 


i = 2,  3.  . . .,  N-1  [Eq  Al] 


However,  internal  shears  can  be  expressed  as  follows  in  terms  of  story 
displacements: 


where  E 
'i 

Wi.Wi^l 


I 

6' 


V2EI. 

(w  -w 

hj  ’ ’ 


) i = 2,3,.  . .,N-1 


-i 

6 


12EIi,^ 

^+1 


= modulus  of  elasticity 
= moment  of  inertia  of  column  at  level  i 
= story  height 

= Story  displacement  at  i^"  and  i^"  + 1 level 


[Eq  A2] 


(w,j_j_1-Wj)  i - 2,3,.  . . ,N-1 


12  = stiffness  parameter  for  the  fixed-end  columns 

6 = number  of  equivalent  columns  (i.e.,  there  are  three 

bays  with  interior  column  stiffness  twice  that  of 
the  exterior  columns). 

To  simplify  the  expressions  and  the  subsequent  derivation,  it 
is  assumed  that  h^  is  equal  to  the  story  height  (h).  Thus, 
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[Eq  A3] 


X I * 

vl  = 6( ^)(w.-w,  J i = 2,3,.  . ..N-1 

1 11-1 


f 12EI.., 

''i+1  " 1 = 2,3,.  . .,N-1 

h 

The  equilibrium  expressions  for  the  frame  then  become: 


Ft  = 6(^^3){-I.w._^  + 

[Eq  A4] 

i = 2,3,. 

. .,N-1 

f{  = v|  - = 6(^){(I^+l2)wi  - I2W2} 

[Eq  A5] 

" ''n  " 

[Eq  A6] 

For  the  shear  wall  system  (Figure  A2),  the  equations  of  equilibrium 
for  the  lateral  force  distribution  can  be  written  as  follows: 


F^  = i = 3,4,5 N-2  [Eq  A7] 

where  F^  = portion  of  design  lateral  force  acting  on  the  shear 
wall  at  floor  i 

= portion  of  shear  in  shear  wall  at  floors  i and  i+1, 
respectively. 

However,  the  internal  shear  on  the  wall  can  be  expressed  in  terms  of 
the  internal  bending  moment  (M^)  of  the  wall  as  follows: 


+ hV^  - = 0 


= 0 


[Eq  A8] 


or 
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[Eq  A9] 


and 


f;  = (-M^_i+2M5-M^+^)/h 


[Eq  AlO] 


From  the  beam-bending  formula  (Eq  All) 


d w _ M 


[Eq  All] 


and  central  difference  considerations,  it  was  determined  that  the  bend- 
ing moment  at  floor  i (M^. ) can  be  expressed  as: 


-El! 


^ ' -;;r<“i+r2«i*“i-i) 


[Eq  A12] 


Thus,  the  lateral  force  on  the  structure  may  be  expressed  in  terms  of 
the  structure  story  displacements  as  follows: 


i; 


5 '7?  “3  - * (6^)*'1»  [E'’  '"3] 

h I,  I,  I, 


"i  =-«-r  'o'  'i 


i: 


l^+i? 


i; 


EjS  jS 

1 'tI  "4  - 204)«3  " - 2('4)«1>[E<I  A"'] 

n I2  I2  I2  I2 


Ed  d 


d 


d ,+d 


‘i  'i 

d d 

- 2('4r)«i-l  ^ «i-2'tEq  A15] 
M M 


i = 3,4,5,.  . . ,N-2 
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•^N-l  ''N-1 


\/5  = — [iii  {_2w  + (4+-Jlb^w  _ 2(1+- -~--)w 

''n  .3  ^ ^ N-2 

''  'n-1  'n-1 


I 

X N-2  , 

,s  V3^ 


I 


[Eq  A16] 


N-1 


El 

F^  = (w  -2w  +w  ) 

''n  ^3  ^'^N  ^Vl  '^N-2^ 


[Eq  A17] 


The  above  set  of  equations  for  and  F^.  (i  = 1,2,.  . .,N)  are  the  equa- 
tions that  govern  the  behavior  of  the  structural  model  due  to  a lateral 

distribution  of  forces  over  the  model.  The  final  equilibrium  expres- 

f s 

sions  for  the  structure  are  obtained  by  adding  F|  and  F^.  for  i =1,2, 

3,.  . .,N. 

Thus 


F.  = f!  I-  F^  = H.(w)  i = 1,2,3 N [Eq  A18] 

where  F^  = total  external  lateral  load 
f s 

F^,F^.  = as  defined  above 

H(w)  = a function  of  the  stiffness  of  the  elements  of  the  struc- 
ture and  the  absolute  lateral  displacements. 

Dead  and  live  load,  which  are  also  assumed  in  the  model,  are  con- 
sidered further  in  the  Design  Considerations  section.  A series  of  runs 
for  typical  building  configurations  was  made  to  examine  the  response 
and  the  moments/forces  in  frame-shear  wall  elements.  The  results  ob- 
tained compared  favorably  with  those  of  Khan  and  Sbarounis,®  whose  work 
provided  the  basis  for  this  check  on  the  structural  model. 

Shear  expressions  for  the  shear  wall  were  included  in  initial 
consideration  of  the  model.  Although  this  unnecessarily  complicated 


® F.  R.  Khan  and  J.  A.  Sbarounis,  "Interaction  of  Shear  Walls  and 
Frames,"  Proceedings  ASCE  Structural  Division^  ST3  (June  1964). 
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the  model,  it  did  provide  some  insight  into  the  relative  magnitude  of 
shear  wall  bending  story  displacement  and  shear  wall  shear  story  dis- 
placement. For  massive  shear  walls  in  low  buildings,  shear  wall  shear 
story  displacement  dominates  shear  wall  bending  story  displacement. 
However,  for  buildings  of  four  or  more  stories  and  even  moderately  heavy 
shear  walls,  shear  wall  bending  becomes  the  dominant  contributor  to 
story  displacement. 


Design  Considerations 

The  computer  program  simulates  a structure's  design  for  a given 
load  requirement  to  provide  consistent  information  for  different  sets 
of  load  requirements.  The  simulated  design  is  generated  to  meet  a 
set  of  criteria.  These  criteria,  which  require  that  the  column  and 
shear  wall  resistances  equal  or  exceed  the  force  demands,  are  as 
follows: 

1.  Columns  shall  be  sized  for  ultimate  strength  considerations 
for  axial  and  bending  moment. 

2.  Columns  shall  be  checked  for  tension  from  overturning  moment. 

3.  Columns  shall  be  checked  for  shear  strength. 

4.  Shear  walls  shall  be  sized  for  elastic  strength  considerations 
in  shear. 

5.  Shear  walls  shall  be  checked  for  tension  from  overturning 
moment. 

Examination  of  the  set  of  forces  on  the  columns  and  shear  wall 
will  facilitate  understanding  of  the  design  simulation  process.  For 
the  frame  (Figure  A3) 

J 6EL^  N N 

■^i  = (w.-w.  ,)  + ^ W.(w.-w.  J ± ^ F.(j-i+l)h(w.-w.  ,)/6. 

h^  1 1 I j=i  J 1 1 I .]  11-11 


f 6EL  N N 

^ W.w,  ± I F.(jh)w^  6/^  i = 2,3,.  . .,N 
‘ h ' j=i  o=i  J ' ' 


[Eq  A19] 


where  the  first  term  is  the  bending  moment  induced  in  a fixed-ended 
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N N 

R=jSW  ±2  F.(j-in)h/8, 

j=i  J j=j  J I 

i = l,2,...,N 


Figure  A3.  Column  axial  load  relationship  for  design  moments. 
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column  from  story  displacement,  the  second  term  is  the  bending  moment 
induced  by  the  dead  and  live  load  acting  eccentrically  on  the  column, 
and  the  third  term  is  the  compression  or  tension  induced  by  the  over- 
turning moment  acting  eccentrically  on  the  column  with  a moment  arm 
of  6,^..  The  expressions  for  shear  in  the  frame  are  as  follows: 


12El! 


,N 


f 12EI 


[Eq  A20] 


For  axial  loads  on  the  frame,  the  following  expression  is  apparent 
from  Figure  A3: 

X N N 

P?  = I W.  ± I F.(j-i+l)h/6.  i = 1,2 N 

j=i  j=i 


[Eq  A21] 


J=i  " J 
For  the  shear  wall  (Figure  A4) 


hN 


-ElS 


N 


2^(w^^1~2w.+w._i)  + SGN  ABS(  J 6Wj(w^-w^^-|) 

N ^ 

i I F.(o-i+l)h(w.-w.  ,)/6.) 

j=i  j 

i = 2,3,.  . 


[Eq  A22] 


El 


,N-1 


M,  = - y ~(w^-2w, ) + SGN  ABS(  T 6W.w, 

^ h^  2 1 j 1 

N 

± I F jhw,/6,), 
j=l  ^ ' ' 


N 


M 


0 


where  again  the  first  term  is  due  to  bending  moment  induced  in  the 
shear  wall  resulting  from  story  displacement,  the  second  term  is  bend- 
ing moment  induced  by  dead  and  live  loads  acting  eccentrically  on  the 
wall,  and  the  third  term  is  the  compression  or  tension  induced  by  the 
overturning  moment  acting  on  the  wall  eccentrically  with  an  arm  c'‘  6^.. 

The  dead  load  and  overturning  moment  increase  the  internal  moment 
on  the  frame  and  shea*  wall  members.  The  signs  of  the  first,  second. 
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and  third  terms  of  Eqs  A19  and  A22  are  in  agreement. 

The  expressions  for  shear  in  the  shear  wall  are  as  follows: 


El- 


I 


I? 

''i  “ -7  <-”1+1  (2-^)  «i  - 0+2-^)  w,_,  + 

n 1^  I, 

[Eq  A23] 


El: 


V=  = (4w,-W2) 


= 


El 


N-1 


tWf^-2W(^_l  + Wj^_2> 


For  axial  loads  on  the  shear  wall,  the  following  expression  is  ap- 
parent from  Figure  A4: 

- N N 

P-  = I 6W.  ± ^ F.(j-i+l)h/6.  i = 1.2 N [Eq  A24] 

j=i  j=i  ' 

The  design  algorithm  currently  considers  three  types  of  struc- 
tural systems  with  three  load  distribution  configurations;  a ductile 
moment-resisting  space  frame  (K  equals  0.67),  a dual  bracing  system 
with  a ductile  moment-resisting  frame  (K  equals  0.80),  and  a box  sys- 
tem (K  equals  1.33).  For  the  system  where  K equals  0.67,  100  percent 
of  the  dead,  live,  and  lateral  loads  on  the  structure  are  carried  by 
the  frame.  In  the  cases  where  K equals  0.80  or  1.33,  the  frame  is 
assumed  to  carry  25  percent  of  the  dead,  live,  and  lateral  loads,  and 
the  shear  wall  system  is  designed  to  resist  100  percent  of  the  dead, 
live,  and  lateral  loads.  It  is  important  to  note  that  in  these  two 
cases,  the  load  on  the  frame  members  is  determined  on  the  basis  of  100 
percent  dead,  live,  and  lateral  load.  The  internal  forces  in  the  mem- 
ber are  then  scaled  to  25  percent,  and  the  resistance  for  an  adequately 
designed  member  must  be  larger  than  this  value.  This  is  a somewhat  dif- 
ferent procedure  from  establishing  the  total  member  internal  forces  for 
25  percent  of  the  applied  loads  and  proceeding  to  design  the  member  on 
this  basis.  In  effect,  the  frame  and  shear  wall  system  are  designed 
together  instead  of  individually  for  certain  percentages  of  load. 
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The  lateral  displacements  of  frame  and  shear  wall  systems  are  compatible, 
resulting  in  a somewhat  different  distribution  of  moment  and  shear  forces 
between  the  frame  and  shear  wall  elements  than  if  an  individual  member 
design  procedure  was  employed.  This  simplied  approach  should  provide 
designs  comparable  to  those  obtained  by  individually  designing  the  frame 
and  shear  wall  members.  Individual  design  of  frame  and  shear  wall  members 
would  have  proved  more  cumbersome  for  this  initial  model  and  certainly 
not  convenient  for  defining  the  structural  configurations  in  terms  of 
ratios  of  moments  of  inertia  of  structural  elements. 

The  frame  design  is  considered  critical  in  the  simulated  design. 

For  even  moderate  ratios  of  shear  wall  to  frame  moments  of  inertia, 
the  shear  wall  system  easily  meets  the  design  requirements. 

Columns  within  the  model  are  sized  on  the  basis  of  ultimate 
strength  considerations  for  columns  under  both  axial  load  and  bending 
moiient.  From  Winter, “ the  ultimate  compressive  and  moment  values  are 
determined  from  the  interaction  diagram  of  Figure  A5. 

For 


Pt*"  " (Ag+A5)y(0.85  f\t) 

where  A^  =•  area  of  tension  steel 

I 

A^  = area  of  compression  steel 

b = column  width 

t = column  depth 

fy  = yield  strength  of  steel 

f^  = yield  strength  of  concrete, 

and 


[Eq  A25] 


6.  Winter  et  al..  Design  of  Concrete  Structures  (McGraw-Hill,  1964). 
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Figure  A5.  Bending  and  axial  load-rectangular  sections. 

Reprinted  with  permission  of  the  American  Concrete 
Institute,  from  C.  S.  Whitney  and  Edward  Cohen, 
"Guide  for  Ultimate  Strength  Design  of  Reinforced 
Concrete,"  Journal  of  the  American  Concrete  Insti- 
tute, Vol  28,  No.  5 (November  1956),  p 455. 
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e 


[Eq  A26] 


where  e = the  eccentricity  ratio 
= ultimate  moment  capacity 
= ultimate  axial  force  capacity  of  the  column, 

the  column  design  algorithm  consists  of  the  following  steps: 

Step  1.  Compute  the  column  eccentricity  from  the  ratio  of  in- 
ternal forces  M and  P: 

e = ” [Eq  A27] 

Step  2.  If  e/t  is  less  than  0.1,  set  e equal  to  0.1. 

Step  3.  Compute  the  column  reduction  factor  R from 

R = 1.07  - 0.008  (h/0. 3b)  [Eq  A28] 

where  h = the  column  height 
b = the  column  width 

Step  4.  Compute  the  point  of  the  intersection  of  the  curves  p^m 

and  e/t  = constant  and  determine  the  interaction  diagram 

column  design  parameters  (Kp  and  from  Figure  A5.  A 

ratio  of  column  effective  depth  to  width  ((})  ) of  0.85  has 

m 

been  selected.  The  ultimate  strength  of  the  column  for 
axial  load  becomes 

P„  • [Eq  A29] 

and  the  ultimate  strength  of  the  column  for  bending  be- 
comes 


M 


of  the  column 


M = P e = K ({)  bt^f' 
u u m^m  c 


[Eq  A30] 
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step  5.  For  the  column  the  following  design  checks  are  made: 

a.  If  is  less  than  M,  the  column  dimensions  are 
increased. 

b.  If  is  greater  than  TM,  the  column  dimensions  are 
decreased.  T is  a convergence  factor  varying  from 
1.5  to  infinity  depending  on  the  extent  to  which  the 
member  resistance  is  allowed  to  exceed  the  require- 
ment. 

c.  If  is  less  than  P,  the  column  dimensions  are  in- 
creased. 

d.  If  P^  is  greater  than  TP,  the  column  dimensions  are 
decreased. 

e.  For  shear,  if  <(>  2/7  b^  is  less  than  V,  the  column 

dimensions  are  increased. 

2 

f.  For  tension,  if  4p  f^b  is  less  than 

N N 

I F.(j-i+l)h/6.  - I W. 
j=i  ^ j=i 

the  column  dimensions  are  increased. 

If  any  of  the  column  dimensions  are  changed  in  Step  5,  the  new 
structural  configuration  is  computed  and  the  process  repeated.  The  pro- 
cess usually  takes  several  iterations  to  converge.  Once  satisfactory 
convergence  has  been  obtained,  the  design  has  been  simulated. 

In  the  shear  wall  design  algorithm,  the  wall  configuration  is  as- 
sumed to  have  fixed  dimensions  and  the  design  parameter  is  the  percent 
steel  in  the  wall.  The  wall  provides  shear  resistance  and  resistance 
to  overturning.  The  nominal  shear  stress  in  the  shear  wall  is  given  by 

Vu  = 4>s(2/f^  + 0.5  p^fy)A^  [Eq  A31] 

where  (j)^  = capacity  reduction  factor  (0.85) 

= shear  resistance 

' 1 

fj,  = compression  strength  of  concrete 


Pg  = percent  steel  in  the  shear  wall 
fy  = yield  strength  of  steel 

= area  of  the  shear  wall. 

The  shear  stress  in  the  wall  cannot  be  greater  than  this  value. 

If  the  value  is  exceeded,  the  percentage  of  steel  in  the  wall  is  in- 
creased. If  the  resistance  is  exceeded  for  a certain  convergence 
option,  the  percentage  steel  in  the  wall  is  decreased. 

The  overturning  moment  is  checked  by  requiring  the  wall  resistance 
to  exceed 

= 0.5  PgfyAV2.0  [Eq  A32] 

where  P^  = overturning  moment  resistance 

P->fy>  and  A^  = as  defined  above. 

5 y 

Here  again,  percentage  steel  in  the  wall  is  increased  for  a design  not 
exceeding  P^.  Based  on  these  expressions,  the  design  is  simulated 
for  a given  structural  configuration  and  lateral  force  requirement. 

The  design  algorithm  may  be  further  complicated  by  constraining  the 
design  frame  and  shear  wall  moments  of  inertia.  It  is  possible, 
and  in  the  past  has  been  useful,  to  require  a frame  member  at  level  i 
to  have  a percentage  of  the  moment  of  inertia  of  a frame  member  at 
level  1.  It  is  also  convenient  to  have  constant  moment  of  inertia 
shear  walls  over  N stories. 

The  degree  to  which  the  simulated  design  agrees  with  the  actual 
design  is  an  indication  of  accuracy  of  the  simulation  model.  However, 
the  designs  may  not  agree  perfectly  because  (1)  the  simulation  design 
is  devoid  of  detail  and  built  around  a very  simple  model  that  behaves 
only  approximately  like  the  real  structure,  and  (2)  the  goal  is  rela- 
tive rather  than  absolute  decision  information.  Since  all  the  simu- 
lated designs  are  generated  in  the  same  way,  the  relative  nature  of  the 
decision  information  is  preserved. 
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The  Costing  Algorithms 


Transforming  structural  behavior  into  information  pertinent  to 
the  decision-making  process  requires  algorithms  that  extract  informa- 
tion from  this  behavior  and  transform  it  into  information  that  can  be 
interpreted  in  terms  of  cost.  The  computer  program  for  the  simulation 
model  contains  three  basic  cost  algorithms— a facility  damage/repair 
cost  algorithm,  a facility  strengthening  cost  algorithm,  and  a facility 
replacement  cost  algorithm.  The  rationale  behind  the  development  of 
each  of  these  algorithms  is  discussed  to  provide  insight  into  the 
generation  of  information  for  decision  making.  Whenever  cost  is 
referred  to,  it  is  assumed  that  the  cost  is  the  construction  cost  for 
the  structure. 

Damage/Repaiv 

The  facility  damage/repair  cost  algorithm  is  based  on  the  computa- 
tion of  a quantity  called  the  damage  factor.  This  damage  factor  relates 
both  the  code  base  shear  requirement  for  the  facility  and  the  struc- 
ture's behavior  to  the  amount  of  damage  that  structural  and  nonstruc- 
tural  components  might  be  expected  to  sustain  if  the  design  code 
requirements  are  exceeded.  The  damage  factor  is  thus  directly  related 
to  the  difference  in  code  requirements. 

It  is  well  known  that  structures  designed  under  different  levels 
of  code  requirements  either  for  different  seismic  zones  or  at  different 
times  experience  damage  when  subjected  to  seismic  activity  exceed- 
ing the  design  requirement,  i.e.,  when  the  demand  on  the  structure 
exceeds  its  capacity.  Some  of  this  experience  is  quantified  in 
Whitman' s“  study  of  the  effects  of  earthquakes  of  different  magnitudes 
on  structures  designed  in  accordance  with  the  Uniform  Building  Code 
(UBC).^^  It  is  also  well  known  that  the  inelastic  behavior  of  the 

R.  V.  Whitman,  Methodology  and  Pilot  Application,  R74-15,  Report 
No.  10,  Structures  Publication  385  (Department  of  Civil  Engineering, 
MIT,  July  1974). 

Uniform  Building  Code  (International  Conference  of  Building  Of- 
ficials, 1973). 
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materials  in  a structure  and  the  action  of  both  structural  and  non- 
structural  elements  provide  the  structure  with  a reserve  energy;  these 
actions  are  often  not  fully  accounted  for  in  the  design  because  of  a 
lack  of  knowledge  of  their  precise  behavior.  Blume“  has  considered 
damage  to  structures  using  a method  termed  the  Reserve  Energy  Technique. 
This  method  relates  the  hysteretic  elastic  energy  capacity  of  a struc- 
ture under  load  to  the  true  combination  elastic-inelastic  energy  capac- 
ity available  in  the  structure.  Using  this  concept,  it  is  possible 
to  relate  the  code  design  requirement  demand,  the  current  facility 
capacity,  and  the  capacity  a structure  may  be  required  to  have  under  a 
higher  demand  to  the  damage  a structure  may  be  expected  to  see  under 
higher  demands.  For  a structure  that  behaves  as  an  elastoplastic,  one- 
dimensional model  with  a ductility  ratio  y. 


A. 

A 


[Eq  A33] 


where  A = a measure  of  capacity  of  a structure  before  collapse 

Ag  = a measure  of  capacity  of  a structure  in  the  elastic  range. 
From  Figure  A6,  equating  the  energy  under  the  demand  "deflection"  curves 
provides  an  expression  for  the  structure's  ductility  in  terms  of  the 
demand  0^  on  the  structure  and  the  capacity  of  the  structure; 


li  = I [(5^)^  1 1] 


[Eq  A3A] 


Once  the  structure's  capacity  has  been  determined  (in  the  case  of 
the  simulation  model,  the  capacity  is  merely  the  design  base  shear)  and 
the  subsequent  demand  on  the  structure  is  identified  (for  the  simulation 
model  this  is  an  anticipated  base  shear  the  structure  may  be  called  upon 


Effects  Pvediotion  Guidelines  for  Structures  Subjected  to  Ground 
Motion,  JAB-99-115,  UC-11  (URS/John  A.  Blume  and  Associates,  Engi- 
neers, July  1975). 
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Elastic  Demand  Model 

Inelastic  Copocity  Model 

Figure  A6.  Demand  and  capacity  energy  models. 
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to  withstand),  Eq  A34  defines  the  extent  the  reserve  energy  will  be 
called  upon  to  withstand  the  demand.  The  ductility  y provides  a rela- 
tive measure  of  the  extension  of  the  structure  behavior  into  the  inelas- 
tic region.  Since  it  is  merely  a relative  measure,  a formula  is  needed 
to  translate  the  measure  into  an  absolute  damage  variable  for  the  struc- 
ture. Blume  suggests  the  rollowing  relationship: 


P ^ 1 

1 ^ y ^ [Eq  A35] 
1 ^ y 


where  = structure  damage  factor 

y = required  ductility  to  meet  the  demand 

u = ultimate  ductility  of  the  structure 

k = an  "economic  factor." 

For  purposes  of  uniformity  in  description,  u^  is  written  as  calibration 
factor  c^  and  k is  written  as  calibration  factor  Cg.  Thus,  DF^j^^  be- 
comes 


DF 


SNS 


0 

y-1 


M&) 


c,-l 
1 


11  s 1 

1 - y 5 c 
S y 


1 


[Eq  A36] 


u is  a relative  measure  of  the  capacity  of  a structure  to  the  point 

r* 

of  collapse.  A value  of  0 implies  that  the  demand  on  the  struc- 
ture was  below  the  design  capacity,  and  for  an  elastic  design,  no  inelas- 
tic reserve  energy  was  required.  For  DF^nj5  values  greater  than  0 but 
less  than  1,  some  inelastic  reserve  energy  being  used  to  resist  the 
demand  resulted  in  damage  to  the  structure.  For  values  of  1» 
the  structure  has  collapsed.  In  terms  of  the  computer  simulation 
model,  Dg  and  are  the  structure's  required  and  designed  base 
shear,  k and  y^  are  calibration  constants,  and  is  an  absolute 

measure  of  structural  damage  to  both  structural  and  nonstructural  com- 
ponents. 

The  damage  factor  DF^j^^  is  not  the  only  possible  measure  of 
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structural  and  nonstructural  damage.  Steinbrugge^**  and  others  have 
indicated  that  there  is  a high  degree  of  correlation  between  relative 
story  displacement  and  nonstructural  damage  (e.g.,  cracks  in  walls, 
glass  breakage,  etc.).  Whereas  provides  an  overview  of  absolute 

structural /nonstructural  damage  for  a one-dimensional  elastoplastic 
model  of  the  structure,  a second  damage  factor  examines  the  structure 
behavior  in  more  detail  to  arrive  at  another  relative  measure  of  damage 
for  nonstructural  elements  of  the  structure.  For  the  simulated  design 
under  the  design  load,  the  relative  story  displacement  is  merely  the 
difference  in  story  displacements  between  the  top  and  bottom  of  the 
story.  For  the  initial  design  subject  to  a higher  demand  requirement 
than  it  was  originally  designed  for,  the  elastic  relative  story  dis- 
placement can  be  computed.  Utilizing  an  expression  similar  to  the  one 
above  for  ductility  gives  the  following  expression  for  the  nonstruc- 
tural damage  factor: 


where  DF^^^  = partial  nonstructural  damage  factor 

A = elastic  relative  story  displacement  for  higher  demand 
= story  drift  criteria 
C2»C4=  calibration  constants. 

Since  both  DF^^^  and  DF^j^^  relate  to  damage  from  seismic  loads,  they  are 
interrelated.  Preselected  weighting  factors  are  used  to  weight  the 
effects  of  the  two  factors  in  obtaining  the  overall  damage  factor  ex- 
pression. The  expression  for  the  total  damage  factor  is  thus 

DF  = w-jOF^j^^  + ^38] 


1 


A < 0 

0 ^ A/D^  < C3 
C3  ^ A/D^ 


[Eq  A37] 


K.  V.  Steinbrugge  et  al . , Studies  in  Seismicity  and  Earthquake 
Dajnage  Statistics,  COM-71-00053  (National  Oceanic  and  Atmospheric 
Administration,  Department  of  Commerce,  1969),  Appendix  A. 
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where  DF 
DF 
DF 


= total  damage  factor 
= overall  design  structure  damage  factor 

NS 

^ weighting  factors. 


relative  story  displacement  damage  factor 


For  purposes  of  this  study,  w^  equals  0.75  and  W2  equals  0.25.  DF  is 
based  upon  four  calibration  constants  that  are  determined  from  real 
data,  and  the  weighting  factors  quantify  the  relationships  between  the 
two  damage  factors. 

The  cost  of  damage/repair  is  proportional  to  DF.  The  structural 
damage/repair  costs  for  the  i^"  floor  (SC^.)  are  written 


SC.  = CSS  . DF3^3 


AR. 


[Eq  A39] 


where  CSS 
AR, 


= Structural  costs  per  square  foot  of  facility 
- area  of  the  i^^’  floor  of  the  facility 


DFsi^S  ' defined  above. 


,th 


The  nonstructural  damage/repair  costs  for  the  i floor  (NSC^. ) are  written 
ten 


NSC.  = CNS 


DF 


AR, 


[Eq  A40] 


where  CNS  = nonstructural  costs  per  square  foot  of  facility 
DF,AR^  = as  defined  above. 

The  total  structural  and  nonstructural  damage/repair  cost  for  a given 
facility  is  merely  the  summation  over  the  floor  levels.  The  rationale 
behind  these  expressions  for  cost  is  based  on  the  difficulty  in  sepa- 


rating the  structural  and  nonstructural  damage  in  DF 


SNS* 


Thus,  both 


SC  and  NSC  contain  indeterminate  nonstructural  and  structural  costs. 

Strengthening 

The  facility  strengthening  model  is  based  on  the  quantity  of  stiff- 
ness that  must  be  added  to  the  facility  to  bring  its  displacements  into 
agreement  with  those  of  the  facility  designed  for  the  higher  require- 
ments. That  is,  the  added  stiffness  (in  matrix  form)  is 
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[Eq  A41] 


where  K^.  = the  stiffness  matrix  for  the  original  structure 

f the  stiffness  matrix  for  a newly  designed  facility  under  the 
hi||ier  requirement. 

Stiffness  is  of  the  form  CEI/Jl^  where  C is  a constant,  E is  the  material 
modulus,  I is  the  section's  moment  of  inertia  (the  only  terntwhich 
varies),  and  Z is  the  length  of  the  dimension.  The  square  root  of 
the  I value  is  proportional  to  area  and  is  taken  proportional  to  the 
cost  of  the  added  material  to  the  facility,  thus  defining  relative 
costs.  In  order  to  "normalize"  the  cost  for  strengthening  over  the 
facility,  the  added  stiffness  is  divided  by  the  average  stiffness  for 
the  original  simulated  facility  design.  This  is  accomplished  by  comput- 
ing  the  average  of  diagonal  terms  of  K^,  and  dividing  into  each  of  the 
diagonal  terms  of  In  a parallel  argument  to  that  used  in  the  dam- 
age/repair model,  a stiffening  factor  for  both  structural  and  nonstruc- 
tural  components  is  defined  in  terms  of  the  following  ratio  of  stiffness 
expressions  for  the  i^  story: 

stiffening  factor  =/ KggCi.i)/!^^  [Eq  A42] 

where  Kggli.i)  = i^”  diagonal  term  of  the  strengthening  stiffness  matrix 
K^.  = average  stiffness  of  the  original  simulated  design 

(average  of  the  diagonal  terms  of  K^). 

Structural  costs  are  proportional  to  the  structural  strengthening  factor 
defined  by 

SFs  = C5/ KggCi.D/K^  [Eq  A43] 

where  SF^  = structural  strengthening  factor 

Cg  = calibration  constant 

/ Kgg(i^!i)/K^.  = as  defined  above. 
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Nonstructural  costs  are  proportional  to  the  nonstructural  strengthening 
factor  defined  by 

SF^S  = Cg/K^5(i,i)/K.  [Eq  A44] 

where  = nonstructural  strengthening  factor 

Cg  = calibration  constant 

/ K,,(i.i)/I<,  = as  defined  above. 

b b I 

The  total  strengthening  factor  (SF)  is  merely  the  sum  of  SFg  and  SFj^jg*. 

SF  = SFg  + SFj^g  [Eq  A45] 

The  cost  of  strengthening  is  proportional  to  the  sum  of  the  indi- 
vidual strengthening  factors  times  the  facility  area  times  the  cost 
per  unit  area  of  the  original  facility.  Thus,  the  structural  strengthen- 

ing  cost  for  the  i'^"  floor  is  written 

SC.  = CSS  • SFg  • AR.  [Eq  A46] 

4*  K 

and  the  nonstructural  strengthening  cost  for  the  i^'*  floor  is  written 

NSC.  = CNS  • SF^3  • AR.  [Eq  A47] 

The  total  structural  and  nonstructural  co't  for  a given  facility  is 
merely  the  summation  over  the  floor  levels. 

Replaoment 

The  facility  replacement  model  is  based  on  measuring  the  increase 
in  material  needed  in  a simulated  design  for  higher  requirements.  Since 
the  lateral-force-resisting  elements  such  as  columns  and  shear  walls 
generally  require  greater  material  quantities  to  withstand  greater 
lateral  forces,  the  increase  in  area  of  column  elements  and  the  increase 
in  steel  ratio  tor  shear  wall  elements  are  taken  to  be  proportional  to 
the  cost  of  facility  replacement. 
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The  facility  replacement  factor  for  both  structural  and  nonstruc- 
tural  elements  is  defined  by 

RF  = C7(b^d^-b^d^)  + Cg(p^-PQ)  [Eq  A48] 

where  RF  = replacement  facility  factor 

Cy,Cg  = calibration  constants  for  frame  and  shear  wall  quantities 
respectively 

b^jd^  = simulated  column  design  width  and  depth,  respectively, 
for  the  original  facility 

b.|,di  = simulated  column  design  width  and  depth,  respectively, 
for  the  replacement  facility 
= simulated  shear  wall  percentage  steel  for  the  original 
facility 

p^  = simulated  shear  wall  percentage  steel  for  the  replace- 
ment facility. 

Available  data  indicate  that  82  percent  of  the  replacement  cost 
can  be  attributed  to  structural  costs  and  the  remaining  18  percent  to 
nonstructural  costs. Thus,  the  incremental  cost  IC  for  the  i story 
may  be  written  as  follows: 

IC  = RF  • AR.  • CSS/0.82  [Eq  A49] 

The  total  structural  cost  then  becomes 

SC.  = 0.82  • IC  + AR.  • CSS  [Eq  A50] 

and  the  total  nonstructural  cost  becomes 

NSC.  = 0.18  • IC  + AR.  • CNS  [Eq  A51] 

R.  V.  Whitman,  Methodology  and  Pilot  Applioation^  R74-15,  Report  No. 
10,  Structures  Publication  385  (Department  of  Civil  Engineering, 

MIT,  July  1974). 
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Again,  the  total  structural  and  nonstructural  cost  for  a given  facil- 
ity is  merely  the  summation  over  the  floor  levels. 

Summary 

Perhaps  the  most  useful  cost  information  is  provided  by  the  aver- 
age of  the  floor  level  costs.  The  computer  program  not  only  provides 
this  information  but  also  the  cost  as  a percentage  of  the  original 
facility  cost,  as  used  in  the  example  in  Chapter  4.  Knowing  the  ori- 
ginal facility  cost  provides  a basis  for  comparing  costs  to  repair, 
strengthen,  or  replace  over  a variety  of  different  facility  configura- 
tions. 
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APPENDIX  B: 


MODEL  ;balibration 


Damage/ Repair  Algorithm 

Four  constants  must  be  determined  For  calibration  of  the  algorithm 
for  damage/ repair- -two  for  DF^j^j.  and  two  for  DF|^g.  In  the  expression 
for  DF^j^2  A36]i,  c-j  represents  the  building's  ul tiiiiate  ducTfl i ty  and’ 
C2  is  an  "economic  factor."  Section  3.3  of  the  American  Concrete  In- 


stitute (ACI)  report  entitled  P.esponae  of  Multistory  Concrete  Struc- 
tural Lateral  Forces^ ^ indicates  that  a recommended  minimum  ductil- 
ity factor  for  reinforced  concrete  buildings  in  earthquake  areas  is 
4.0  to  6.0.  For  estimating  damage,  then,  using  an  ultimate  ductility 
factor  of  4.0  would  be  conservative.  It  should  be  noted,  however,  that 
determination  of  a building  system's  ultimate  ductility  is  difficult, 
and  any  number  assigned  to  broad  classes  of  buildings  may  be  subject  to 
considerable  error. 

An  MIT*^  study  of  four  specially  designed  reinforced  concrete  build- 
ings to  determine  a maximum  ductility  factor  for  buildings  designed 
for  several  UBC  requirements  and  a "super  zone"  (twice  UBC  Zone  3) 
subject  to  a 0.27  g ground  acceleration,  obtained  lower  values  for 
maximum  ductility  in  most  cases.  Table  B1  presents  a portion  of  the  MIT 
data.  These  data  seem  to  indicate,  for  example,  that  a six-story  con- 
crete moment-resisting  frame  building  designed  for  UBC  3 would  have  a 
ductility  factor  of  3,8  when  subjected  to  a 0.27  g peak  ground  acceler- 
ation. Moving  from  UBC  0 design  to  UBC  2 design,  in  general,  decreases 
the  maximum  ductility  inherent  in  the  building.  The  last  two  columns 


Response  of  Multistory  Concrete  Structures  to  Lateral  Forces,  ACI 
SP-36  (American  Concrete  Institute,  1973). 

R.  V.  Whitman,  Methodology  and  Pilot  Application,  R74-15,  Report  No. 
10,  Structures  Publication  385  (Department  of  Civil  Engineering, 

MIT,  July  1974). 
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Table  B1 


Maximum  Ductility  Factors 
With  0.27  g Peak 


Caused  by  Ground  Motion 
Acceleration* 


Maximum  Ductility  Factor 


Building** 

Stories 

UBC  0 

UBC  2 

UBC  3 

S 

'^max 

0 

max 

r." 

r “T 

CMRF 

6 

2.8 

4.3 

0.5 

CMRF 

11 

DIs] 

bji 

2.7 

1.2 

3.2 

0.5 

CSW 

11 

1876  i 

im 

2.2 

1.7 

5.8 

3.3 

CSW 

17 

0*2 1 

r3"i*« 
li -i 

2.0 

1.0 

- 

- 

* From  R.  V.  Whitman,  Methodology  and  Pilot  Application^  R74-15, 
Report  No.  10,  Structures  Publication  385  (Department  of  Civil 
tr.ji nearing,  mIT,  July  ly74). 


**  CSW 

= Concrete  Shear  Wall. 

:mfr 

= Concrete  Moment-Resisting  Frame. 

' □ 

= probable  partial  or  total 

collapse. 

1 — 1 
c _ 1 

= possible  partial  or  total 

collapse. 
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of  Table  B1  give  the  mean  and  standard  deviations  of  the  maximum  duc- 
tility ratios  of  the  concrete  moment- resisting  frame  buildings  of  six 
and  11  stories  and  the  concrete  shear  wall  building  of  11  stories  (con- 
figurations in  Table  B1  corresponding  to  configurations  developed  in 
the  model)  over  the  UBC  0,  UBC  2,  and  UBC  3 design  strategies.  Since 
will  be  determined  over  the  number  of  building  stories  without  re- 
ference to  a moment-resisting  frame  or  a shear  wall  configuration,  to 
be  consistent,  c^  must  also  be  determined  on  this  basis.  For  the  six- 
story  concrete  moment- resisting  frame,  the  mean  and  standard  deviation 
over  the  first  three  design  strategies  are  4.3  and  0.5,  respectively. 
For  the  11 -story  concrete  moment-resisting  and  shear  wall  configura- 
tions, the  mean  and  standard  deviation  of  the  maximum  ductility  ratio 
were  calculated  from  the  values  in  the  table  to  be  4.5  and  2.5,  respec- 
tively. If  ultimate  ductility  is  assumed  to  correlate  with  "probable 
partial  or  total  collapse"  of  Table  B1 , then  both  4.3  and  4.5  are  low, 
since  both  are  weighted  toward  "possible  partial  or  total  collapse"-- 
the  former  more  so  than  the  latter.  In  addition,  the  data  to  be  used 
in  determining  C2  are  probably  weighted  toward  the  lower  design  strat- 
egy requirements  (UBC  0,  UBC  2),  also  implying  an  increased  mean  maxi- 
mum ductility  factor.  For  simplicity,  a value  of  ultimate  ductility 
of  5.0  has  been  selected  for  all  building  configurations  of  Table  B1 . 
The  uncertainties  in  the  data  do  not  allow  a more  refined  estimate. 

A preliminary  determination  of  C2  was  made  from  the  data  on  the 
San  Fernando  earthquake  compiled  at  MIT,*®  which  provides  perhaps  the 
most  extensive  data  base  available  on  damage  to  structures  from  seismic 
activity.  From  Blume's  work,*®  a structure's  ductility  is  expressed  in 
terms  of  its  design  capacity  and  the  demand  placed  on  it  (Eq  A34).  The 
damage  factor  can  be  expressed  in  terms  of  the  ratio  of  the  demand  on 


^*®  R.  V.  Whitman,  Damage  Probability  Matrices  for  Prototype  Buildings, 
R73-57,  Report  No.  8,  Structures  Publication  380  (Department  of 
Civil  Engineering,  MIT,  October  1973). 

* ® Effects  Prediction  Guidelines  for  Structures  Subjected  to  Ground 
Motion,  JAB-99-115,  UC-11  (URS/John  A.  Blume  and  Associates,  Engi- 
neers, July  1975). 
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the  structure  to  the  structure's  design  capacity.  Here,  p is  assumed 
to  represent  the  structure  ductility  at  demand  for  a given  capacity 
Cg.  The  expression  for  the  damage  factor  can  be  written  as  follows: 


log  (DF3^3)  = C2  log  (^)  1 < p < c^  [Eg  Bl] 


If  c^  and  p are  known  in  terms  of  the  demand  to  capacity  ratio,  then  a 
relationship  between  the  damage  factor  for  various  demands  and  C2  can 
be  developed.  If  the  damage  factor  can  be  correlated  with  a demand,  the 
value  of  C2  can  be  determined. 

The  MIT  data  base  on  the  San  Fernando  earthquake  is  expressed  in 
terms  of  damage  probability  matrices.  Damage  probability  matrices 
(Figure  Bl ) express  the  damage  state  of  a structure  in  terms  of  the 
percentage  of  structures  in  that  state  (P^gj)  for  a particular  building 
configuration  and  a particular  demand  level  expressed  in  terms  of  MMI. 
The  damage  state  can  be  correlated  with  the  damage  ratio  central  value, 
as  in  Table  B2.  The  MMI  can  be  correlated  with  the  nominal  peak  ground 
acceleration,  as  in  Table  B3.  Thus,  the  damage  a structural  configura- 
tion may  be  expected  to  experience  (defined  by  the  percentage  of  struc- 
tures in  a particular  damage  state)  can  be  related  to  a value  of  nominal 
peak  ground  acceleration  (a  reasonable  measure  for  defining  seismic 
activity).  A particular  damage  probability  matrix  can  be  summarized  at 
a given  intensity  by  computing  a mean  damage  ratio  (MDR).  This  weighted 
average  over  the  damage  states  is  defined  as  follows: 

n. 


_ 1 


•I 


mDR 


I 


1 

■I  i=l 


DRii 


[Eq  B2] 


where  nj  = total  number  of  buildings  in  a particular  category  sub- 
ject to  ground  motion  intensity  I 

DR.r  = damage  ratio  (cost  of  repair  to  cost  of  replacement)  for 
' ^ th 

the  i building  in  a particular  category  subject  to 
ground  motion  intensity  I. 

This  number  provides  a measure  of  the  ratio  of  cost  to  repair  to  cost 
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Figure  B1 . Form  of  damage  probability  matrix. 
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Table  B2 


Correlation  of  Damage  State  to  Damage  Ratio  Central  Value 


Damage  Damage  Ratio* 

State  Central  Value,% 


0 

1 

2 

3 

4 

5 

6 

7 

8 


0 

0.1 

0.5 

2.0 

5.0 

10.0 

30.0 

100.0 

100.0 


*Damage  Ratio  = 


cost  of  repair 
cost  of  replacement 


Table  B3 

Correlation  of  MMI  to  Nominal  Peak  Ground  Acceleration 


MM  I 


Nominal  Peak  Ground  Acceleration,  g 


IV 

V 

VI 

VI. 5 
VII 
VIII 
IX 
X 


0.007 

0.015 

0.030 

0.050 

0.090 

0.200 

0.500 

1.020 


7fi 


4 


of  replacement  for  the  particular  building  configurations  that  are  used 
to  develop  the  matrix.  For  computing  the  mean  damage  ratio  is 
taken  to  be  the  damage  factor  DF5|^|5* 

The  capacity  of  the  structure  is  assumed  to  be  1968  SEAOC  Zone  3 
design.  Then,  the  base  shear  for  the  design  (V)  may  be  expressed 

V = KCZW  [Eq  B3] 

where  K = structural  system  type  coefficient 

C = structural  system  dynamic  coefficient 
Z = seismic  zone  factor 
W = total  weight  of  structure. 

Its  capacity  (u)  can  be  expressed  as 

U = 1.4V  = 1.4KCZW  [Eq  B4] 

To  use  the  information  from  the  MIT  study,  it  is  necessary  to  cor- 
relate response  spectra  with  MMI.  For  3 percent  damping,  Figure  B2 
correlates  the  elastic  response  spectra  bounds  with  MMI  based  on  the 
correlation  of  MMI  with  nominal  peak  ground  acceleratin'  orovided  in 
Table  B3.  A response  spectrum  approach  is  used  to  establish  the  demand 
on  the  structure  for  given  peak  nominal  ground  acceleration  values. 

For  reference  purposes,  the  1968  SEAOC  Zone  3 capacity  and  the  1.0  g 
response  spectrum  demand  values  are  provided  in  Table  B4  over  the  four-, 
seven-,  and  ten-story  configurations,  two  dimensional  ratios,  and 
three  structural  systems.  The  capacity  of  the  facility  can  be  read 
from  the  table  and  the  demand  scaled  from  the  value  corresponding  to  a 
1.0  g response  spectrum. 

Table  B5  presents  the  required  structure  ductility  based  on  the 
1968  SEAOC  Zone  3 capacity  and  the  response  spectrum  demand  determined 
for  MMIs  of  VI,  VII,  and  VIII.  The  variation  in  correlation  of  MMI 
with  nominal  peak  ground  acceleration,  the  selection  of  3 percent  damp- 
ing for  spectra  amplification  factor  computation,  and  the  accuracy  in 
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Correlation  of  elastic  response  spectra  bounds  with  MMI.  SI  conversion 
factor:  1 in. /sec  = 25.4  mm/sec. 


Table  B4 


Capacity  Demand  for  Various  Structural  Configurations 

Four  Story 


K 

T 

Dimensional 

Ratio 

Capacity 
Cs=  1.4,KCZW 

Demand  * 
Ds  = CSgW 

0.67 

.4 

1.0 

70,656 

1,032,300 

0.67 

.4 

1.6 

113,049 

1 ,651 ,680 

0.80 

.843 

1.0 

71,130 

1,008,000 

0.80 

.667 

1.6 

123,081 

1,612,800 

1.33 

.843 

1.0 

127,122 

825,600 

1.33 

.667 

1.6 

219,966 

1,320,960 

Seven  Story 


0.67 

.7 

1.0 

101,415 

1,722,240 

0.67 

.7 

1.6 

162,267 

2,755,584 

0.80 

1.476 

1.0 

103,293 

■ 1,686,300 

0.80 

1.167 

1.6 

178,737 

2,698,080 

1.33 

1.476 

1.0 

186,444 

1,383,960 

1.33 

1.167 

1.6 

322,617 

2,214,336 

Ten  Story 

0.67 

1.0 

1.0 

128,037 

2,416,050 

0.67 

1.0 

1.6 

204,858 

3,865,680 

0.80 

2.108 

1.0 

131,022 

2,370,000 

0.80 

1.667 

1.6 

226,716 

3,792,000 

1.33 

2.108 

1.0 

237,426 

1,929,300 

1.33 

1.667 

1.6 

410,835 

3,086,880 

* 


In  the  response  spectrum  approach,  the  ratio  of  C to  the  response 
spectrum  acceleration  (C/S  ) is  defined  as  follows: 

® Stories 


k 


Four 


Seven 

Ten 

0.897 

0.885 

0.803 

0.790 

0.607 

0.590 

0.67 

0.80 

1.33 


0.9.' 

0.840 

0.640 
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Table  B5 


Model  Ductility  Based  on  Capacity  and  Demand 


Freq 

MMI 

mm! 

Mf^I 

y 

Hz 

VI 

VII 

VIII 

VI 

VII 

VIII 

Four  Story 

2.50 

70,656 

89,810 

269,430 

598,734 

1.31 

7.77 

36.40 

2.50 

113,049 

143,696 

431 ,088 

957,974 

1.31 

7.77 

36.40 

1.19 

71,130 

67,536 

201 ,600 

443,520 

0.95 

4.52 

19.94 

1.50 

123,081 

140,314 

420,941 

935,424 

1.15 

6.35 

29.38 

1.19 

127,122 

55,315 

165,120 

363,264 

0.59 

1.34 

4.58 

1.50 

219,966 

114,924 

344,771 

766,157 

0.64 

1.73 

6.57 

Seven  Story 

1.43 

101,415 

139,501 

421 ,949 

930,010 

1.45 

9.16 

42.55 

1.43 

162,267 

223,202 

675,118 

1,488,015 

1.45 

9.16 

42.55 

0.68 

103,293 

65,766 

202,356 

455,301 

.70 

3.11 

10.21 

0.86 

178,737 

129,508 

404,712 

863,386 

.76 

3.06 

12.17 

0.68 

186,444 

53,974 

166,075 

373,669 

.54 

1.04 

2.51 

0.86 

322,617 

106,288 

332,150 

708,588 

.55 

1.03 

2.91 

Ten  Story 

1.0 

128,037 

137,715 

422,809 

918,099 

1.08 

5.95 

26.21 

1.0 

204,858 

220,344 

676,494 

1,468,958 

1.08 

5.95 

26.21 

0.47 

131,022 

63,990 

189,600 

450,300 

.62 

1.55 

6.41 

0.60 

226,716 

128,928 

398,160 

872,160 

.66 

2.04 

7.90 

0.47 

237,426 

52,091 

154,344 

366,567 

.52 

.71 

1.69 

0.60 

410,835 

104,954 

324,122 

709,982 

.53 

.81 

1.99 

an 


reading  values  from  the  elastic  bound  response  spectrum  plots  may  re- 
sult in  as  much  as  25  percent  error  in  some  of  the  computed  demand 
values  in  the  table.  This  is  within  the  accuracy  of  the  model  for 
damage/repair. 

The  which  is  taken  equivalent  to  the  mean  damage  ratio  of 

the  MIT  study,  is  presented  in  Table  B6  for  five  to  seven  story  and 
eight  to  13  story  buildings  for  MMIs  of  VI,  VII,  and  VIII.  Table  B7 
presents  the  overall  damage  probability  matrices  from  the  MIT  study. 

The  structures  making  up  the  data  sample  were  constructed  prior  to  1933 
or  after  1947.  The  post-1947  data  are  presented  along  with  the  combina- 
tion pre-1933  and  post-1947  data.  In  Table  B8  the  calibration  constant 
^2  is  presented  for  the  four-,  seven-,  and  ten-story  buildings  over  the 
building  configuration  and  MMIs  of  VI,  VII,  and  VIII.  Blanks  occur  in 
the  table  when  the  computed  ductility  is  either  greater  than  5.0  or 
less  than  1.0.  Certain  cases  examined  showed  high  ductilities  but 
very  low  damage  ratios.  This  is  unreasonable  and  biases  the  values  of 
to  the  high  side.  Values  cf  C£  less  than  5 are  the  only  ones  con- 
sidered. 

Based  on  averages  over  each  story  configuration  assuming  C2  to  be 
independent  of  MMI,  the  following  values  of  Cg  and  their  standard  de- 
viations are  selected: 


Four  Story 

C2  = 2.66 

(^c2 

= 0.68) 

Seven  Story 

C2  = 2.78 

Q 

0 

ro 

= 1.67) 

Ten  Story 

C2  = 2.18 

^^c2 

= 0.85) 

For  the  portion  of  the  damage/repair  algorithm  related  to  rela- 
tive story  displacement  (Eq  A37),  c^  is  taken  to  be  1.0,  i.e.,  nonstruc- 
tural  damage  related  to  the  relative  story  displacement  is  a linear 
function  of  the  story  displacement.  There  are  no  data  and/or  guidance 
available  that  indicate  otherwise. 

To  complete  calibration  of  this  portion  of  the  algorithm,  c^  is 
determined  from  some  sparse  data  relating  relative  story  displacement 
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Table  B6 


Mean  Damage  ratio  (DFcwc)  for  Concrete  Buildings  of  Five  to  Seven 
and  Eight  to  13  Stories  for  MMIs  of  VI,  VII,  and  VIII 


Five  to  Seven  Story 


MMI 

Post-1947 

★ 

All 

VI 

.0002 

.0002 

VII 

.0105 

.0268 

VIII 

.0267 

.0267 

Eight  to  13  Story 

MMI 

Post- 1947 

•^''SNS 

All 

VI 

.0015 

.0011 

VII 

.0043 

.0162 

VIII 

.0963 

.0963 

* Includes  all  da*-a  available  in  the  study  regardless  of  con 
sfruction  era. 


82 


' B5I.  AVAILABLE  COPY 

Table  B7 


Summary  of  Damage  Matrices,  Replacement  Cost  Version* 


* R.  V.  Whitman,  et  al..  Damage  Statistics  for  High-Rise  Buildings 
in  the  Vicinity  of  the  San  Fernando  Eartl^uake,  R73-24,  Report  No. 
7,  Structures  Publication  363  (Department  of  Civil  Engineering, 
MIT,  April  1973). 

**  Data  for  damage  states  expressed  in  percentage'.  MDRs  and  standard 
deviations  are  not  expressed  in  percentage, 
t Mean  damage  ratio. 
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Table  B8 


Cp  Computed  on  the  Basis  of  MIT  San  Fernando 
Earthquake  Data  and  DF^j^g 


Post-1947 

All 

yi 

yi_ 

3.33 

3.33 

Four 

3.33 

3.33 

Story 

- 

- 

2.59 

2.59 

Post-1947 

All 

Vi 

VI 

3.90 

3.90 

Seven 

3.90 

3.90 

Story 

• 

** 

•* 

Post-1947 

All 

VI 

yi 

1.66 

1.74 

Ten 

1.66 

1.74 

Story 

- 

- 

Post-1947 

All 

Post-1947 

VII 

VII 

VIII 

. 

• 

• 

- 

35.64 

28.31 

- 

1.85 

1.47 

32.66 

2.68 

2.13 

- 

Post-1947 

All 

Post-1947 

VII 

VII 

VIII 

4.40 

3^49 

- 

6.87 

5.45 

- 

- 

- 

3.72 

.93 

.74 

4.90 

Post-1947 

All 

Post-1947 

VII 

VII 

VIII 

. 

•m 

- 

- 

- 

2.75 

2,08 

- 

4.05 

3.06 

- 

- 

- 

1.33 

- 

~ 

1.68 
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to  a mean  damage  ratio  from  the  San  Fernando  earthquake.  Figure  B3 
presents  these  data.  It  should  be  noted  here  that  the  data  in  Figure 
B3  include  both  structural  and  nonstructural  damage.  However,  for 
purposes  of  this  model,  the  structural  damage  is  assumed  to  be  insig- 
nificant so  that  the  damage  ratio  is  very  nearly  proportional  and 
on  the  same  order  of  magnitude  as  DFj^^*  In  reference  to  Figure  B3, 
the  MIT  study  states,  "the  results  appear  to  fall  into  two  groups: 

One  group  with  larger  damage  from  two  buildings  with  many  stiff  and 
brittle  partitions  not  isolated  from  the  structura''  frame,  and  a 
second  group  with  smaller  damage  from  buildings  in  which  there  were 
either  few  partitions  or  with  flexible  and/or  isolated  partitions."^® 
The  data  are  weighted  heavily  toward  nonstructural  damage.  By  fitting 
a straight  line  to  the  data,  DFj^^  may  be  expressed  in  terms  of  the 
relative  story  displacement  as  follows: 


DFj^3  = 0.06486  + 0.0016 


[Eq  B5] 


where  6 = relative  story  displacement  in  feet. 
Thus,  solving  for  6,  realizing  that 


D = 0.06h 


[Eq  B6] 


where  h is  the  story  height  in  feet,  and  substituting  into  Eq  A37 
yields 


^3  = 


DFj^g  - 0.0016 

T:0648)DF^ 


Table  B9  tabulates  c^  for  various  values  of  DF^^g 


[Eq  B7] 
for  a building  with 


story  height  of  10.0  ft  (3.0  m).  From  the  values  in  this  table. 


R.  V.  Whitman,  Damage  Probability  Matrices  for  Prototype  Buildings ^ 
R73-57,  Report  No.  8,  Structures  Publication  380  (Department  of 
Civil  Engineering,  MIT,  October  1973). 
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earthquake.  From  R.  V.  Whitman,  Damage  Probabili-ty  Matinees  for  Prototype 
Buildings,^  R73-57,  Report  No.  8,  Structures  Publication  380  (Department  of 
Civil  Engineering,  MIT,  October  1973).  SI  conversion  factor:  1 ft  = 0.3048  m. 


Table  39 


Correlation  of  DF|^^  With 


<^3 

0.0100 

21.60 

0.0500 

24.90 

0.1000 

25.31 

1.0000 

25.68 
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[Eq  B8] 


c-  = 24.37,  0^  = 1.88,  and  COV  (c^)  = 0.08. 
IS  taken  as  or 


= 24.37. 

This  completes  calibration  of  the  damage/repair  algorithm. 

Strengthening  Algorithm 

For  the  strengthening  algorithm,  calibration  requires  the  deter- 
mination of  two  constants— Cg  and  Cg— that  translate  model  behavior 
into  strengthening  costs  (Eqs  A43  and  A44).  The  ratio  of  nonstructural 
to  structural  strengthening  can  be  taken  as  the  average  of  the  data 
on  the  four  facilities  shown  in  Table  BIO. 

For  Letterman,  Hays,  Oakland  Naval,  and  Charleston  Naval  Hospi- 
tals several  strengthening  schemes  with  different  costs  were  proposed. 
Table  Bll  provides  the  data  based  on  crude  estimates  of  costs  for  the 
four  hospital  configurations  with  one  workable  strengthening  scheme 
for  each.  The  approximate  ratio  of  the  capacity  of  the  strengthened 
facility  to  the  original  facility  is  computed  along  with  a ratio  of 
strengthening  cost  to  initial  facility  costs.  The  figures  are  very 
approximate  based  on  the  best  information  that  could  be  derived  from 
the  reports  associated  with  each  of  these  facility  strengthening  pro- 
posals.^* The  capacity  of  the  hospitals  is  often  different  in  two  di- 
rections. The  capacity  is  taken  to  be  the  minimum  in  this  study.  The 
facility  capacity  ratios  do  not  correlate  well  with  the  cost  ratios, 
due  in  part  to  the  peculiarity  of  each  strengthening  scheme. 


^ Seismic  Safety  Investigation  of  Charleston  Naval  Hospital  at 
Charleston,  South  Carolina  (John  A.  Blume  and  Associates,  Novem- 
ber 1974);  Seismic  Safety  Investigation  of  Oaklayid  Naval  Hospital 
at  Oakland,  California  (John  A.  Blume  and  Associates,  November 
1974);  P.  J.  Richter,  Concepts  for  Strengthening  Letterman  and 
Hays  Army  Hospitals — Task  4C,  Aseismic  Design  Criteria  of  Existhig 
Hospital  Facilities,  Memorandum  No.  7338-4C-3241  (Agbabian  As- 
sociates, June  1974). 
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Table  BIO 


Ratio  of  Nonstructural  Costs  to  Structural  Costs  for 
Strengthening  of  Four  Critical  Facilities 


Hospital 


CNS/CS 

0.154 


Letterman 
Oakland  Naval 
Charleston  Naval 


0.067 

0.024 
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strengthening  Cost  Data 


Facility 

Cost  to  -6 

Cost 

Capacity 

Hospital 

Strengthen  (10  ) 

Ratio 

Ratio 

Letterman 

11.9 

.74 

6.5 

Hays 

4.6 

.29 

5.9 

Oakland 

10.5 

.64 

3.6 

Charleston 

.8 

.06 

2.8 

To  arrive  at  calibration  constants  for  the  model,  two  sets  of 

data  from  Table  Bll  are  considered.  Letterman,  Hays,  and  Oakland  are 

considered  as  one  case  for  arriving  at  strengthening  costs  for  the 

four-  and  seven-story  models,  and  Letterman  and  Oakland  as  a second 

* 

case  for  arriving  at  strengthening  costs  for  ten-story  models.  This 
division  was  used  because  the  Letterman  and  Oakland  Hospitals  corre- 
late with  ten- story  models,  whereas  Hays  Hospital  is  better  represented 
by  a seven-story  model  in  the  current  study.  However,  it  would  be 
dangerous  to  conclude  on  the  basis  of  the  sparse  data  in  Table  Bll  that 
low-level  buildings  have  strengthening  costs  significantly  less  than 
those  for  high-rise  buildings  (particularly  since  Hays  Hospital  is  only 
two  to  three  stories  shorter  than  Letterman  or  Oakland).  The  simula- 
tion model  does,  however,  seem  to  indicate  a trend  in  this  direction; 
thus,  the  data  for  low-level  structures  provide  a weighting  of  strength- 
ening costs  on  all  models.  A weighted  average  is  therefore  the  most 
acceptable  alternative  for  the  low-level  structures  at  this  time. 

The  cost  to  strengthen  the  ten-story  facility  model  from  the 
simulation  is  taken  to  be  69  percent  with  a 5.1  facility  capacity 
ratio  (an  average  of  the  values  in  Table  Bll  for  Letterman  and  Oakland 
Hospitals).  The  cost  to  strengthen  the  four-  and  seven-story  facility 
model  from  the  simulation  is  taken  to  be  56  percent  with  a 5.3  ■facil- 
ity capacity  ratio. 

To  calibrate  the  model  for  strengthening  for  ten-story  configura- 
tions, all  ten-story  cost  data  for  which  k equals  1.33  are  averaged 
for  arbitrary  Cg  and  Cg.  The  constants  Cg  and  Cg  are  scaled  to  the 
above  percentage  strengthening  cost  for  the  facility  capacity  ratio. 

For  four-  and  seven-story  configurations,  all  seven-  and  ten-story  cost 
data  for  which  k equals  1.33  are  averaged  for  arbitrary  c^  and  c,  with 

3 D 


* Charleston  Naval  Hospital  was  not  considered  further  in  this  study 
because  its  comparatively  low  facility  capacity  ratio  indicates 
a special  case  of  strengthening  not  compatible  with  the  other 
three  cases.  In  fact,  significant  structural  strengthening  was 
not  recommended  ir  the  case  of  the  Charleston  Naval  Hospital. 
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Cg  and  Cg  scaled  as  in  the  above  case.  Table  B12  provides  the  strength 
ening  coefficients. 

Replacement  Algorithm 

The  replacement  algorithm  calibration  requires  the  determination 
of  constants  c^  and  Cg  in  Eq  A43.  Data  for  calibration  of  the  model 
are  taken  from  two  sources.  Appendix  B of  a report  of  the  SEAOC  Ad 
Hoc  Committee  on  Costs  of  Design  for  Earthquakes^^  states  that  the  cost 
for  design  and  construction  of  a reinforced  concrete  structure  for  UBC 
Zone  2 is  102  percent  of  the  cost  for  design  and  construction  for  UBC 
Zone  0.  For  UBC  Zone  3,  this  publication  gives  the  cost  increase  as 
5 percent.  A second  source  of  data  computes  cost  figures  from  the  de- 
sign of  a six-story  and  an  11-story  reinforced  concrete  apartment 
building  for  UBC  Zones  0,  1,  2,  3,  and  S (Zone  S requirements  are  twice 
UBC  Zone  3 requirements).^^ 

A linear  extrapolation  of  these  data  for  a critical  facility  for 
1.0  g response  spectrum  acceleration  using  the  SEAOC  and  MIT  data 
yields  an  average  increase  in  replacement  costs  of  74.8  and  30.3  per- 
cent, respectively.  The  74.8  percent  increase  is  probably  excessive, 
particularly  in  shear  wall  buildings  where  the  configuration  may  remain 
basically  the  same  at  low  and  high  levels,  but  the  percentage  steel 
in  the  shear  wall  is  increased.  For  this  study,  the  four-  and  seven- 
story  configurations  are  based  on  data  for  the  six-story  concrete 
building  in  the  MIT  study,  and  the  ten-story  configuration  is  based  on 
the  data  for  the  11-story  concrete  building. 

For  four-story  configurations,  extrapolation  of  the  MIT  data 
based  on  a Zone  3 requirement  and  capacity  and  demand  values  at  SEAOC 
Zone  3 and  a 1.0  g response  spectrum  indicated  replacement  costs  to  be 

Percent  Inavease  in  Design  Inspection  and  Construction  Cost^ 

Report  of  the  Ad  Hoc  Committee  on  Costs  of  Design  for  Earthquakes 
(SEAOC). 

R.  V.  Whitman,  Methodology  and  Pilot  Application,  R74-15,  Report  No 
10,  Structures  Publication  385  (Department  of  Civil  Enqineei  i>^q, 
MIT,  July  1974). 
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Table  B12 

Simulation  Model  Coefficients 


C2 

^7 

^8 


Four  Story 

Seven  Story 

Ten  Story 

5.00 

5.00 

5.00 

2.66 

3.38 

2.18 

24.37 

24.37 

24.37 

1.00 

1.00 

1.00 

1.812 

1.812 

1.78 

0.0317 

0.0317 

0.0312 

0.00998 

0.00349 

0.00155 

76.109 

26.650 

11.801 
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on  the  order  of  124  percent.  A similar  consideration  of  seven-  and 
ten-story  facilities  provided  replacement  costs  of  136  percent  and  131 
percent,  respectively.  All  three  structural  configurations  (K  equals 
0.67,  0.80,  and  1.33)  were  considered  in  computing  these  costs,  c^ 
and  Cg  were  determined  by  equating  the  percentage  increase  in  costs  to 
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The  ratio  of  Cg  to  c^  was  determined  by  requiring  the  wall  reinforce- 
ment steel  quantity  to  be  equal  to  the  frame  reinforcement  steel  qual- 
ity. This  provided  realistic  percentages  of  shear  wall  steel  rein- 
forcement. Table  B12  presents  the  values  for  c^  and  Cg. 

It  should  also  be  noted  here  that  structural  costs  (CSS)  are  taken 
to  be  82  percent  of  the  total  cost  increase  (TCIN).  The  data  of 
Figure  B4  provide  some  support  for  this  percentage.  Thus  the  struc- 
tural cost  increase  is  .82  TCIN  and  the  nonstructural  cost  increase 
is  .18  TCIN,  where 


TCIN  = [c^lb^d^-b^d^)  + Cg(p^-p^)]/0.82  [Eq  BIO] 

The  model  is  calibrated  for  replacement  data. 

A final  computation  /or  the  costs  for  a critical  facility  should  be 
presented.  Five  case  studies  of  Veterans  Administration  hospital  con- 
struction^'* give  the  ratio  of  structural  to  nonstructural  costs  as 
25.2  percent.  The  total  nonstructural  cost  is  given  by 

CNS  = 0.80  TC  = 0.80  (41.03)  = $32.82/sq  ft 

($353.27/m^)  [Eq  Bli] 
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Feasibility  Study — V.A.  Hospital  Building  System^  Research  Study 
Report,  Project  Number  99-R003  (Building  Systems  Development  and 
Stone,  Marraccini  and  Patterson,  October  1968). 


CO 

H 

Z 

UJ 

z 

O 

a. 

2 

8 

< 

oc 

r> 

8 

o: 

H* 

CO 

§ 

Z 

8 

UJ 

3 

O 

Z 

2 

Z 

a 

h“ 

CO 

O 

U 


S5 


CMRF  - Concrete  Moment  Resisting  Frome 
SMF,  SMRF  - Steel  Moment  Resisting  Frame 


Figure  B4.  Relative  cost  premium  for  structural  and 

nonstructural  components.  From  R.  V.  Whitman, 
Methodology  and  Pilot  Applioation,  R74-15, 
Report  No.  10,  Structures  Publication  385  (De- 
partment of  Civil  Engineering,  MIT,  July  1974). 
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and  the  total  structural  cost  by 


CS  = 0.20  TC  = 0.20  (41.03)  = $8.21/sq  ft  ($88.37/m^)  [Eq  B12] 
The  costs  are  averaged  over  the  five  hospitals. 

Sumnary 

The  calibration  data  presented  above  for  the  three  algorithms  can 
be  summarized  as  follows: 

1.  Damage/repair  data  consider  only  story  height  breakdown  (five 
to  seven  stories  and  eight  to  thirteen  stories)  and  building  material 
(steel  or  concrete).  Building  system  and  ratios  of  relative  siffnesses 
are  not  a part  of  the  classification.  For  nonstructural  damage/repair 
data,  there  is  no  breakdown  at  all. 

2.  Strengthening  data  considers  building  system  and  story  height; 
however,  the  available  data  are  for  a K of  1.33  only. 

3.  Replacement  data  provide  no  breakdown  according  to  building 
configuration.  Extrapolation  of  the  data  for  various  story  height 
configurations  results  in  different  calibration  constants  because  of 
the  differing  ratios  of  structure  demand  to  capacity. 

The  model  contains  more  detail  than  can  be  provided  by  the  avail- 
able data.  The  calibration  process  is  kept  simple  so  that  any  uni- 
formity in  variation  of  the  costs  can  be  identified  over  the  wall  to 
frame  stiffness  parameter.  Observations  such  as  general  decrease  in 
strengthening  costs  with  increase  of  wall  to  frame  stiffness  are  a prod- 
uct of  this  'ffort.  As  a result  of  this  averaging  of  the  calibration 
constants  over  the  wall  to  frame  stiffness  configuration,  the  cost 
information  will  be  averaged  over  the  wall  to  frame  stiffness  config- 
uration also.  Thus,  calibration  to  131  percent  replacement  cost  for  a 
ten-story  building  will  show  a variation  from  138  percent  replacement 
cost  for  a frame  structure  down  to  a 109  percent  replacement  cost  for 
a light  shear  wall  structure.  This  variation  over  structural  config- 
uration leaves  the  model  well  within  the  accuracy  of  the  available  data. 
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SYMBOLS 


AR^.  - area  of  the  floor. 

- area  of  shear  wal 1 . 

Ag  - area  of  column  tension  steel. 

I 

Ag  - area  of  column  compression  steel, 

b - column  width. 

c^(i=lj  8)  - model  calibration  constants. 

C - structural  system  dynamic  coefficient. 

I I 

CpC.  - cost  coefficients. 

Cg  - structure  capacity  in  terms  of  base  shear  force. 

CNS  - total  nonstructural  construction  costs  per  square 

foot. 

COV(x)  - coefficient  of  variation  of  x. 

(C/S^)  - response  spectrum  structural  system  dynamic  coef- 

a 

ficient. 

CSS  - total  structural  construction  costs  per  square  foot, 

d^.  - column  effective  depth. 

- story  drift  design  requirement. 

Dg  - demand  on  the  structure  in  terms  of  yield  base  shear 

force. 

dj^.dj^  - structure  story  displacement  at  level  x (at  level  N). 

OF  - total  damage  factor. 

H 
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e 

E 


I 


F,  F. 
h 

H(w) 


I 


c 


XL 

- damage  factor  for  nonstructural  damage  (for  the 
story). 

- damage  factor  for  structural  and  nonstructural 
damage. 

- damage  ratio  (cost  of  repair/cost  of  replacement) 

XL 

for  the  structure  in  a particular  category  sub- 
ject to  ground  motion  intensity  I. 

- column  eccentricity  ratio, 

- modulus  of  elasticity. 

- an  event. 

- yield  strength  of  concrete. 

- yield  strength  of  steel. 

- portion  of  design  lateral  force  acting  on  the  frame. 

- total  external  lateral  load. 

- portion  of  design  lateral  force  acting  on  the  shear 
wall . 

- design  lateral  force  at  height  x. 

- structure  overturning  moment  forces. 

- story  height 

- structural  height  from  ground  to  x. 

- a function  of  the  stiffness  of  the  elements  of  the 
structure  and  the  absolute  lateral  displacements. 

moment  of  inertia  of  structure  model  exterior  column. 


I 


f 


I 


s 


IC 

k 

K 


Kss(iJ) 


"^DRI 

M 


NSC. 

P 


- moment  of  inertia  of  structure  model  frame. 

- moment  of  inertia  of  structure  model  shear  wall. 

- incremental  cost. 

- an  economic  factor  (Blume  damage  model). 

- structural  system  type  coefficient. 

- existing  structure  stiffness  matrix. 

- average  stiffness  of  the  original  simulated  design. 

- interaction  diagram  column  design  parameters. 

- strengthened  structure  stiffness  matrix. 

- strengthening  parameter  stiffness  matrix. 

- i^^  diagonal  term  uf  the  strengthening  parameter 
stiffness  matrix. 

- length  of  dimension. 

- mean  damage  ratio  at  MMI  I. 

- bending  moment. 

- moment  frame  at  level  i. 

- moment  in  shear  wall  at  level  i. 

- ultimate  moment  capacity  of  a column. 

- total  number  of  buildings  in  a particular  category 
subject  to  ground  motion  intensity  I. 

- nonstructural  damage/ repair  costs  for  the  i level. 

- percentage  of  steel  in  one  column  face. 
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percent  steel  in  shear  wall 


OSI 


P(E) 

P 

u 

R 

RF 


SNS 


SCi 


SF 

SF. 


NS 


t 

T 

TC 

TCIN 

U 

V 


'b 

,f 


- column  compression  and  bending  design  parameter. 

- percentage  of  structures  in  a particular  damage 
state. 

- external  force  on  plane  at  i^^  level. 

- external  force  on  shear  wall  at  i^^  level. 

- probability  of  occurrence  of  event  E. 

- ultimate  axial  force  capacity  of  a column. 

- column  reduction  factor. 

- replacement  factor  for  structural  and  nonstructural 
elements. 

- spectrum  acceleration  level. 

'fh 

- structural  damage/repair  costs  for  the  i^^"  level. 

- strengthening  factor  for  nonstructural  elements. 

- strengthening  factor  for  structural  elements. 

- effective  column  depth. 

- convergence  factor. 

- total  structure  construction  cost. 

- total  cost  increase. 

- structure  capacity. 

- structure  design  base  shear. 

- structure  shear  at  base, 
portion  of  shear  in  frame. 

- portion  of  shear  in  shear  wall. 
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- ultimate  column  shear  resistance. 

- structure  yield  base  shear. 

th 

- story  displacement  at  the  i level. 

- weighting  factors  for  determination  of  DF 

- total  weight  of  structure  including  dead  load  and 
live  load  (weight  of  i story). 

- seismic  zone  factor. 

- relative  story  displacement  (for  the  i story). 

- relative  measure  of  structure  displacement  before 
collapse  (single  degree  of  freedom  model). 

- relative  measure  of  structure  displacement  at  elas- 
tic limit  (single  degree  of  freedom  model). 

- relative  story  displacement  in  feet. 

- dimension  factor  for  overturning  moment. 

- ultimate  story  strain. 

- structure  ductility. 

- ultimate  structure  ductility. 

- standard  deviation  of  x. 

- capacity  reduction  factors. 
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